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ABSTRACT 


An  experimental  investigation  was  conducted  to  study  convective  heat  transfer  in 
straight  and  curved  rectangular  channels  of  high  aspect  ratio  that  approximate  plates 
of  infinite  extent.  Experiments  were  performed  at  steady  state  in  the  turbulent  flow 
regime  with  one  wall  held  at  a  constant  heat  flux  and  the  opposite  wall  essentially 
adiabatic.  The  effect  of  curvature  induced  secondary  flow  on  heat  transfer  on  the 
concave  and  convex  w’alls  was  observed  by  comparing  Nusselt  numbers  for  four 
different  configurations  at  several  different  Reynolds  numbers.  Significant  heat  transfer 
enhancement  was  observed  on  the  concave  wall.  Correlations  for  Nusselt  number  as  a 
function  of  Reynolds  number  were  calculated  for  the  cases  studied. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Lord  Rayleigh  [Ref.  1]  was  seemingly  the  first  to  investigate  the  mechanics  of 
flows  with  curved  streamlines.  His  analysis  showed  that  inviscid  curved  streamline 
flows  are  stable  providing  that  circulation  increases  with  radial  distance  (for  flow 
between  concentric  cylinders  this  condition  is  met  when  the  outer  cylinder  spins  more 
quickly  than  the  inner  cylinder). 

Later,  Taylor  [Ref  2J  studied  viscous  flow  between  concentric  rotating  cylinders 
and  showed  that  when  the  inner  cylinder  is  rotated  while  holding  the  outer  cylinder 
stationary,  the  onset  of  instability  is  predicted  by  the  dimensionless  parameter: 

Ld  nr 

Ta  (Tavlor  number)  =  —  J —  £41.3.  (eqnl.l) 

v  ri 

where  U  is  the  tangential  velocity  of  the  inner  cylinder,  d  is  the  gap  between  cylinders, 
r.  is  the  radius  of  curvature  of  the  inner  cylinder,  and  v  is  kinematic  viscosity.  The 
instability  is  manifest  in  the  formation  of  regularly  spaced,  counter-rotating  vortices 
whose  axes  of  rotation  lie  in  the  circumferential  (main  flow)  direction  (see  Figure  1.1). 
Furthermore,  Taylor  concluded  that  these  vortices  form  in  the  laminar  flow  regime  but 
do  not  induce  turbulent  main  flow  [Ref.  3J. 

The  flow  of  a  viscous  fluid  in  a  curved  channel  was  studied  by  Dean  [Ref.  4|.  lie 
developed  an  analytical  solution  for  flow  in  a  curved  channel  approximated  by 
concentric  cylinders  where  the  radius  of  the  inner  cylinder  was  much  greater  than  the 
gap  between  cylinders.  Instead  of  rotating  the  inner  cylinder  Dean  imposed  a  pressure 
gradient  along  the  channel  wall  and  predicted  that  flow  instability  would  occur  when  a 
critical  value  of  a  dimensionless  parameter  similiar  to  the  Taylor  number  was  exceeded: 

Id  nr 

De  (Dean  number)  =  —  V  —  £36,  <eqn  I  .) 

v  ri 

where  d  represents  the  channel  half-width. 
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Between  Rotating  Cylinders 


Gortler  [Ref.  5]  went  on  to  investigate  viscous  flow  in  a  curved  channel  and 
concluded  that  the  (low  instability  predicted  by  Dean  would  take  a  form  analagous  to 
the  Taylor  vortices  of  concentric  cylinder  flow  (see  Figure  1.2).  Regularly  spaced 
counter  rotating  vortices  form  along  axes  of  rotation  parallel  to  the  main  flow 
direction.  Yih  and  Sangster  [Ref.  6]  concluded  that  the  centrifugal  force  coupled  with 
the  velocity  gradient  near  the  concave  wall  caused  instability  due  to  a  stratification  in 
specific  weight  that  led  to  vortex  formation. 

The  secondary  flows  predicted  by  Taylor  and  Gortler  (Taylor-Gortler  vortices) 
have  been  confirmed  experimentally  using  various  flow  visualization  techniques,  and 
hot  wire  anemometry  [Ref.  7).  The  presence  of  these  vortices  in  the  laminar  flow 
regime  led  to  speculation  of  the  elfect  of  this  secondary  flow  on  heat  transfer  over 
curved  surfaces.  Extensive  research  has  been  done  which  yields  good  experimental 
agreement  with  analytical  solutions  for  laminar  flow  in  curved  channels  of  rectangular 
cross  section  [Refs.  8,9,10).  and  for  turbulent  flow  in  a  curved  rectangular  channel  of 
low  aspect  ratio  t ratio  of  the  spanwise  dimension  to  the  channel  height) 
[Refs.  11,12,13).  Research  in  these  areas  indicates  that  there  is  enhancement  of  heat 
transfer  in  these  flow  regimes  when  Taylor-Gortler  vortices  form  near  the  concave  wall. 
However,  information  relating  to  the  fluid  mechanics  and  heat  transfer  in  the  transition 
and  turbulent  flow’  regimes  for  channels  of  high  aspect  ratio  is  very  limited. 

B.  APPLICATIONS 

Studying  heat  transfer  over  a  curved  surface  is  just  one  aspect  of  the  complex 
fluid  dynamics  and  heat  transfer  problems  that  exist  in  an  operating  gas  turbine.  The 
high  temperatures  characteristic  of  a  combustion  engine  constrain  gas  turbine 
performance  because  of  metalurgical  limitations  of  turbine  blades  which  in  turn 
restricts  the  efficiency  of  the  turbine.  Information  gathered  in  this  type  of 
experimentation  could  lead  to  a  better  understanding  of  the  problem  and  facilitate  the 
design  of  gas  turbines  of  higher  efficiency  and  power  to  weight  ratio. 

C.  INTENT  OF  THIS  STUDY 

The  purpose  of  this  investigation  was  to  study  the  effect  of  curvature  on  heat 
transfer  in  a  curved  rectangular  channel  of  large  aspect  ratio  by  gathering  extensive 


14 


data  relating  Nusselt  number  to  Reynolds  number  for  flow  over  walls  of  constant  ueat 
flux  of  straight,  convex  and  concave  geometry.  This  information  will  provide  a  firm 
data  base  for  validating  computer  models  of  the  complex  flow  and  heat  transfer  in 
curved  rectangular  channels  similar  in  geometry. 

Data  for  hydraulically  fully  developed  flow  in  a  straight  rectangular  channel 
heated  either  from  above  or  below  was  taken  to  provide  a  basis  for  comparison  with 
the  flow  over  a  heated  convex  or  concave  wall.  Empirical  correlations  similiar  to  the 
familiar  Dittus-Boelter  correlation  for  Nusselt  number  as  a  function  of  Reynolds 
number  were  calculated  from  the  data  so  that  the  heat  transfer  enhancement  could  be 
described  quantitatively  and  compared  to  previous  data. 

A  secondary  goal  was  to  determine  the  relative  effect  of  buoyancy  on  heat 
transfer  in  a  curved  channel.  The  effect  of  buoyancy  on  heat  transfer  in  forced  flow 
must  be  considered  if  the  effects  of  free  and  forced  convection  are  comparable  because 
the  combined  heat  transfer  is  then  dependent  on  the  direction  of  free  convection  in 
comparison  to  that  of  forced  convection. 
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II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


A  detailed  description  of  the  materials  and  construction  methods  used  to  build 
the  experimental  channel  and  associated  instrumentation  is  contained  in  G.  Galyo's 
thesis  [Ref.  14:  pp.  25-43].  Some  minor  changes  had  to  be  made  to  the  existing 
hardware  and  computer  software  and  are  described  herein. 

A.  DESCRIPTION  Of  THE  APPARATUS 

The  channel  was  built  of  plexiglass,  having  a  channel  width  of  25.4  cm  and  a 
height  of  0.635  cm,  yielding  an  aspect  ratio  of  40:1.  The  channel  was  rigidly  supported 
on  a  wooden  frame  and  insulated  over  its  entire  length  with  three  layers  of  half  inch 
foam  insulation.  The  curved  section  subtends  ISO5  of  an.,  and  the  convex  wall  (inner 
curved  wall)  had  a  radius  of  curvature  of  29.7  cm  (see  Figure  2.U  xir  was  pulled 
through  the  channel  by  a  compressor  joined  to  the  duct  by  two  inch  PVC  pipe  through 
a  throttle  valve. 

Four  30.48  cm  x  25.4  cm  wire  wound  rubber  heaters  laminated  to  copper  plates 
were  mounted  on  the  inside  walls  of  the  channel  to  accomplish  heating  of  the  air.  Two 
heaters  were  mounted  in  the  straight  test  section,  one  on  the  upper  wall  and  one  on  the 
lower,  at  a  sufficient  distance  from  the  channel  entrance  to  insure  that  (low  was 
hydrodynamically  fully  developed  prior  to  flowing  over  the  straight  heated  sections. 
The  two  remaining  heaters  were  installed  in  the  curved  test  section,  one  on  the  concave 
wall  and  one  on  the  convex  wall,  with  the  leading  edges  of  the  heaters  beginning  at  an 
angle  of  101.5°  from  the  curved  section  entrance;  this  allowed  vortices  to  form  prior  to 
flowing  over  the  curved  heated  sections.  The  heaters  were  constant  power  density  type 
capable  of  600  watts  each,  and  could  be  energized  individually  or  in  any  combination 
of  pairs  by  use  of  two  power  supplies. 

Thi  tv  copper-constantan  thermocouples  were  imbedded  in  each  of  the  copper 
plates  at  the  copper  plate-rubber  heater  interface  so  that  local  temperatures  could  be 
measured  (see  Figure  2.2).  Four  thermocouples  were  mounted  laterally  across  the 
channel  entrance  to  measure  inlet  temperature  (see  Figure  2.1).  Similiarly,  four 
thermocouples  had  been  mounted  at  the  e\a  of  the  channel;  however,  experimentation 


Straight  Test  Sections 


Figure  2.1  Experimental  Rectangular  Channel. 
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Compressor 


showed  that  there  was  insufficient  mixing  of  the  air  up  to  the  exit  which  resulted  in 
erratic  outlet  temperature  measurements  [Ref.  14:  p.  67].  Instead,  eighth  inch 
honeycombing  was  installed  at  the  exit  of  the  channel  to  promote  mixing,  and  the  four 
exit  thermocouples  were  placed  in  the  PVC  pipe  near  the  end  of  the  channel  outlet 
transition  bell  (see  Figure  2.1).  Two  layers  of  half  inch  foam  insulation  were  applied  to 
the  PVC  to  a  distance  of  approximately  fourty-eight  centimeters  past  these 
thermocouples  to  minimize  heat  loss  to  the  surroundings  prior  to  measuring  outlet 
temperature.  One  thermocouple  was  placed  downstream  of  the  flow  measuring  orifice 
and  used  to  calculate  the  density  of  the  air  flowing  through  the  orifice.  Three  groups 
of  five  thermocouples  each  were  connected  in  parallel  between  the  layers  of  insulation 
surrounding  the  heated  sections  so  that  losses  to  the  surroundings  through  the 
insulation  could  be  estimated. 

Any  combination  of  eighty  voltages  from  the  thermocouples  and  power  supplies 
could  be  read  by  the  Hewlett  Packard  3054A  automatic  data  acquisition  system  used 
for  this  research.  Coupled  to  the  data  acquisition  system  were  a  Hewlett  Packard  9826 
computer  and  a  printer  which  provided  control  of  the  data  acquisition  system  and 
facilitated  recording,  reducing,  and  printing  of  data. 

A  How  measuring  orifice  was  installed  in  accordance  with  the  ASME  Power  Test 
Code  (Ref.  15],  and  water  and  mercury  manometers  were  used  to  measure  the  pressure 
upstream  of  the  orifice  and  differential  pressure  across  the  orifice  so  that  mass  flowrate 
could  be  calculated. 

B.  EXPERIMENTAL  PROCEDURE 

The  procedures  followed  for  testing  were  the  same  irrespective  of  the  specific  test 
configuration.  The  procedures  are: 

1.  Zero  the  manometers,  shut  the  system  throttle  valve,  then  start  the  compressor. 

2.  Use  the  throttle  valve  to  adjust  the  differential  pressure  across  the  flow 
measuring  orifice  to  that  corresponding  to  the  desired  flowrate  or  Reynolds 
number.  Determining  the  approximate  differential  pressure  corresponding  to  a 
specific  Reynolds  number  takes  some  experience  and  the  use  of  the 
approximation: 
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3.  Energize  the  heater  and  adjust  the  voltage  for  the  approximate  value  of 
(Tom  -  Tj  )  desired.  In  these  experiments  a  value  of  six  to  twelve  degrees 
centigrade  was  used  for  the  temperature  differential,  the  objective  being  to 
maximize  (T  -Tta)  while  limiting  the  highest  local  plate  temperature  to  a 
maximum  of  about  100°C.  Maximizing  this  temperature  differential  lessens  the 
uncertainty  of  the  calculated  Nusselt  number. 

4.  Experience  showed  that  the  channel  reached  steady  state  about  four  hours  after 
first  energizing  the  heater,  and  approximately  one  hour  to  steady  out  between 
consecutive  data  runs  at  different  flowrates.  When  all  of  the  following  criteria 
were  met,  steady  state  was  considered  to  have  been  achieved: 

a.  Chanee  in  the  quantity  (T  -T.  )  of  less  than  two  percent  over  a  ten 
minute  period. 

b.  Change  in  the  heated  wall  temperature  of  less  than  one  percent  over  a  ten 
minute  period. 

c.  Change  in  heated  wall  insulation  temperatures  of  less  than  one  percent  over 
a  ten  minute  period. 

5.  When  the  above  criteria  for  steady  state  are  met,  run  the  data  gathering  and 
reducing  program  (see  Appendix  A)  several  times  in  succession  to  verify  that 
steady  state  has  been  reached.  Three  to  six  consecutive  data  collection  runs  at 
steady  state  meeting  the  criteria  of  (4)  above  are  required  to  determine  an 
average  Nusselt  number  versus  Reynolds  number  relationship.  Generally  the 
run-to-run  variation  in  Nusselt  number  and  Reynolds  number  were  much  less  at 
lower  flowrates. 

6.  When  a  sufficient  number  of  consistent  data  runs  have  been  recorded,  throttle 
to  adjust  flowrate  to  the  next  desired  Reynolds  number  and  adjust  the  heater 
voltage  to  achieve  the  desired  plate  temperature. 

The  atmospheric  pressure,  orifice  ditferential  pressure  and  pressure  upstream  of 
the  orifice  had  to  be  read  from  a  barometer  and  two  manometers  and  entered  into  the 
data  reduction  program  interactively;  these  parameters  were  not  read  by  the  automatic 
data  acquisition  system. 


During  this  study  it  was  noted  that  small  (one  degree)  ambient  temperature 
perturbations  had  a  noticeable  effect  on  temperatures  measured  throughout  the  test 
system;  therefore,  close  control  of  ambient  temperature  oscillations  was  required  so 
that  temperature  changes  in  the  test  system  could  be  minimized  and  consistent  data 
measured  from  run  to  run. 

The  initial  configuration  of  the  test  channel  was  such  that  buoyancy  induced  flow 
opposed  forced  flow  in  the  curved  test  section.  After  data  was  taken  over  the  full 
range  of  flow  that  the  compressor  was  capable  of,  the  rig  was  inverted  and  data  was 
taken  with  buoyancy  induced  flow  aiding  forced  flow.  The  purpose  of  this  aspect  of 
the  study  was  to  determine  whether  natural  convection  contributed  significantly  to  the 
total  heat  transfer  in  the  curved  test  sections. 
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III.  DATA  REDUCTION 


Each  heater  plate  provided  a  constant  surface  heat  flux.  The  heater  plates  were 
energized  by  regulated  D.C.  power  supplies  through  stable  precision  resistors  of  low 
resistance  value  (see  Figure  3.1).  The  voltage  drop  across  both  the  heater  element  and 
the  precision  resistor  were  read  and  recorded  by  the  data  acquisition  system.  Then  the 
total  power  supplied  by  the  heater  was  determined  from  the  relation: 

V.  x  v 

Qp  *  - 1 (eqn  3.1) 

Pr 

where  Rpr  was  a  known  value  of  resistance  included  as  a  constant  in  the  data  reduction 
program.  The  insulated,  unheated  walls  were  treated  as  adiabatic  because  the  heat  loss 
through  them  was  negligible  compared  to  the  total  heat  input  from  the  heater  plate. 
Therefore,  the  bulk  temperature  varied  linearly  from  the  leading  to  the  trailing  edge  of 
the  heated  plate. 

Thermocouple  voltages  were  converted  to  centigrade  temperature  using  individual 
thermocouple  calibration  curves.  All  temperatures  were  recorded  then  combined  to 
calculate  specific  average  temperatures  desired. 

The  actual  heat  loss  to  the  surroundings  was  difficult  to  account  for.  Typically 
the  estimated  total  heat  loss  by  conduction  through  the  insulation  and  by  radiation 
was  only  about  ten  percent  of  the  actual  total  heat  loss.  The  major  difficulty  in 
making  an  accurate  prediction  of  the  heat  loss  was  due  to  not  being  able  to  determine 
the  area  of  the  channel  over  which  heat  was  lost  to  the  surroundings.  For  the  purpose 
of  estimating  the  heat  loss,  the  area  of  the  heated  plate  was  used  as  the  effective  area 
of  loss;  however,  the  plexiglass  section  of  the  channel  adjacent  to  the  heater  plates  acts 
as  extended  heat  transfer  surfaces,  and  the  other  copper  heater  surfaces  in  the  channel, 
which  have  a  low  resistance  to  thermal  conduction,  act  like  heat  sinks.  Therefore, 
instead  of  assuming  that  the  energy  transferred  to  the  air  was  equal  to  the  energy  input 
from  the  plate  minus  the  estimated  energy  lost  to  the  surroundings  over  the  area  of  the 
plate,  the  energy  actually  transferred  to  the  working  fluid  was  calculated  as  follows: 
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Qair  =  mcpAT  (eqn  3.2) 

where  AT  =  Tout  -  Tin 

The  corresponding  plate  average  heat  transfer  coefficient  and  Nusselt  number  based  on 
the  hydraulic  diameter  of  the  duct  were  then  calculated  as  follows: 


h  = 


lair 


APl  <Twh  -  Tb> 
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(eqn  3.3) 


(eqn  3.4) 


where  DhJ.  the  hydraulic  diameter  of  the  channel,  is  a  geometric  constant, 
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4  *  Cross  Sectional  Area  of  the  Duct 
Wetted  Perimeter  of  the  Duct 


(eqn  3.5) 


Data  on  atmospheric  pressure,  orifice  differential  pressure,  and  orifice  upstream 
pressure  was  used  with  the  flow  measuring  orifice  characteristics  and  air  density  based 
on  air  temperature  at  the  orifice  to  calculate  mass  flowrate,  m.  Reynolds  number 
based  on  hydraulic  diameter  was  then  found  by  the  relation: 
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(eqn  3.6) 


where  Hair  was  the  viscosity  of  the  working  fluid  at  the  average  bulk  temperature 
measured  in  the  test  channel.  Similiarly.  Reynolds  numbers  based  on  the  pipe  diameter 
and  channel  height  were  also  determined.  For  the  concave  and  convex  sections  only, 
the  dimensionless  parameter.  Dean  number,  was  calculated  from: 
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(eqn  3.7) 
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IV.  EXPERIMENTAL  RESULTS 
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Several  sets  of  data  were  taken  over  the  approximate  range  of  Reynolds  numbers 
from  6000  to  20,000.  A  computer  program  (see  Appendix  A)  was  used  to  reduce  the 
raw  data  and  calculate  the  quantities  of  interest,  particularly  plate  average  Nusselt 
numbers  for  corresponding  Reynolds  numbers,  both  based  on  the  hydraulic  diameter  of 
the  channel.  Three  to  six  separate  data  sets  at  each  Reynolds  number  were  taken,  then 
averaged  to  give  the  value  of  the  quantities  of  interest. 

After  extensive  data  had  been  taken  for  the  channel  configured  such  that  free 
convection  opposed  the  direction  of  forced  convection,  the  channel  was  inverted  and 
several  sets  of  data  were  taken  for  the  convex  and  concave  plates  with  the  directions  of 
free  and  forced  convection  coincident  so  that  the  effect  of  buoyant  flow  could  be 
determined. 

The  significant  parameters  and  corresponding  uncertainty  for  the  different 
configurations  follows  in  Tables  1  through  4.  Graphic  comparisons  of  plate  averaged 
Nusselt  number  versus  Reynolds  number  are  shown  in  Figures  4.1  through  4.7.  These 
graphs  also  include  calculated  uncertainty  ranges  for  the  data,  the  best  fit  correlations 
for  each  configuration,  data  for  comparison  from  other  sources,  and  some  of  the 
correlations  that  have  been  documented  in  research  similiar  to  this  experimentation. 

A.  HEATING  THE  STRAIGHT  UPPER  PLATE 

A  total  of  sixty-nine  data  pairs  were  taken  over  an  approximate  range  of 
Reynolds  numbers  of  6800  to  20.000  (see  Table  1).  One  means  used  to  analyze  the 
data  was  the  statistics  software.  MINI  TAB  [Ref.  16].  It  was  observed  that  the  data  at 
Reynolds  numbers  below  7000  exhibited  a  relatively  large  deviation  from  the  regression 
in  comparison  to  the  data  for  Reynolds  numbers  above  7000.  This  result  agrees  with 
previous  research  done  by  J.  Wilson  [Ref.  1~:  p.  5<>|  which  revealed  that  transition  How 
in  the  channel  persisted  up  to  a  Reynolds  number  of  approximately  7100;  therefore, 
data  at  Reynolds  numbers  below  7<hm)  was  not  included  in  the  regression  calculated  for 
the  upper  straight  plate. 
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Using  a  least  squares  regression,  a  correlation  for  the  heat  transfer  from  the 
upper  plate  was  determined  to  be: 


Nuhd  *  «024ReM1)-7s 


(eqn  4. 1 ) 


where  Reynolds  and  Nusselt  numbers  are  based  on  the  hydraulic  diameter.  The 
statistical  correlation  of  the  regression  was  0.998.  The  experimental  uncertainty  of  the 
data  was  approximately  ±5%  for  Nusselt  number  and  ±6%  for  Reynolds  number. 

For  comparison,  consider  the  well  known  Dittus-Boelter  correlation  based  on 
hvdraulic  diameter  and  a  Prandtl  number  of  0.71: 


NUhd  =  0.02Reh 


(eqn  4.2) 


At  first  glance  this  small  difference  between  the  correlations  seems  insignificant; 
however,  a  look  at  Figure  4.1  shows  that  the  Dittus-Boelter  correlation  overestimates 
the  actual  Nusselt  number  by  about  thirty  to  thirty-seven  percent  over  its  range  of 
applicability  (Re£  10.000). 

Also  shown  for  comparison  on  Figure  4.1  are  interpolated  values  from  the  Kays 
and  Leung  [Ref.  18]  analytical  solution  for  infinite  parallel  plates.  These  values  too  lie 
outside  the  uncertainty  limits  of  the  actual  data,  but  only  overestimate  Nusselt  number 
by  ten  to  fifteen  percent  when  compared  to  the  finite  fiat  plates  of  this  research. 

Data  previously  taken  by  G.  Galyo  [Ref  14:  p.  5S]  is  included  in  this  plot.  His 
data  agrees  with  that  of  this  research  to  within  the  bounds  of  uncertainty.  It  must  be 
noted  that  the  data  taken  by  Galyo  at  Reynolds  numbers  above  about  17.000  is  less 
certain  than  the  other  data  due  to  the  difficulty  he  encountered  in  recording  an 
accurate  channel  outlet  temperature  at  high  flowrates  [Ref.  14:  p.  67], 

B.  HEATING  THE  STRAIGHT  LOWER  PLATE 

A  total  of  sixty-four  data  pairs  were  recorded  over  an  approximate  range  of 
Reynolds  numbers  of  6000  to  19000  (sec  Table  2).  As  in  the  presious  case  for  the 
upper  plate,  the  data  recorded  for  Reynolds  numbers  less  than  about  "ixm)  were  not 
included  in  the  least  squares  regression  because  of  their  relatively  high  deviation  from 
the  regression.  The  best  fit  regression  correlation  was: 


The  statistical  correlation  of  the  regression  was  0.998,  and  the  experimental  uncertainty 
of  the  data  was  essentially  the  same  as  that  for  the  upper  flat  plate. 

Figure  4.2  shows  that  over  the  full  range  of  testing  there  is  essentially  no 
difference  in  the  results  for  the  two  flat  plate  cases.  Data  previously  taken  by  G. 
Galyo  (Ref.  14:  p.  59]  is  shown  as  well  for  comparison,  and  is  well  within  the 
experimental  uncertainty  of  the  data  recorded  in  this  study  except  at  the  highest 
Reynolds  number  where,  as  mentioned  previously,  G.  Galyo  s  data  was  highly 
uncertain. 

C.  HEATING  THE  OUTER  CURV  ED  (CONCAVE)  PLATE 

Seventy-two  data  pairs  were  taken  over  an  approximate  range  of  Reynolds 
numbers  of  60UU  to  20.0UO  (see  Table  3).  A  least  squares  regression  of  the  data 
resulted  in  the  relation: 

Huj  *  005Rehd°  M  <e9n  4  J> 

The  statistical  correlation  of  the  regression  was  0.999,  and  the  experimental  uncertainty 
was  the  same  as  that  for  the  two  previous  cases. 

Figure  4.3  revealed  that  the  previous  data  reported  by  G.  Galyo  [Ref.  14:  p.  60] 
coincides  well  with  the  data  recorded  in  this  study  except,  again,  at  the  highest 
Reynolds  number  where  G.  Galyo  had  difficulty  determining  duct  outlet  temperature 
accurately.  Also  shown  for  comparison  are  extrapolated  correlations  from  Kreith 
[Ref.  12:  p  1253],  and  Brinich  and  Graham  [Ref.  II:  p.  35]. 

In  contrast  to  the  results  for  the  Hat  plates  there  was  no  significant  deviation 
from  the  regression  even  at  the  lowest  Reynolds  number  tested  for  the  concave  plate. 
This  supports  previous  research  by  J.  W  ilson  [Ref.  1':  p.  5<>|  which  indicated  that  the 
limiting  Reynolds  number  for  transition  How  in  curved  sections  was  less  than  that  for 
flat  plates  and  approximately  equal  to  53<h». 


The  data  for  the  flat  plates  was  used  as  a  baseline  for  comparison  with  the 
concave  plate,  and  heat  transfer  enhancement  for  the  concave  section  was  reflected  by 
an  increase  in  Nusselt  number  of  about  twenty-one  percent  at  a  Reynolds  number  of 
10,000  to  fourteen  percent  at  a  Reynolds  number  of  20,000  (see  Figure  4.4). 

There  is  a  tendency  to  extrapolate  this  data  to  where  the  correlations  for  the 
concave  plate  and  the  flat  plates  converge  in  order  to  estimate  a  Reynolds  number 
where  the  effect  of  curvature  is  no  longer  significant.  For  these  correlations 
convergence  occurs  at  a  Reynolds  number  of  about  205,000  and  a  Nusselt  number  of 
approximately  232;  however,  research  by  Brinich  and  Graham  indicated  that  the 
enhancement  in  heat  transfer  in  a  concave  section  was  significant  even  at  Reynolds 
numbers  as  high  as  250,000. 

D.  HEATING  THE  INNER  CURVED  (CONVEX)  PLATE 

A  total  of  seventy-three  data  pairs  were  recorded  over  a  range  of  Reynolds 
numbers  of  6000  to  20,000  (see  Table  4).  A  least  squares  fit  of  the  data  gave  the 
relation: 

~\d  -  0.017Rehd°'78  (eqn  4.5) 

at  a  statistical  correlation  of  0.997,  and  an  experimental  uncertainty  the  same  as  for  the 
previous  cases. 

Figure  4.5  revealed  that,  as  was  the  case  for  the  concave  plate,  there  was  no 
significant  deviation  from  the  regression  even  at  the  lowest  Reynolds  number  tested.  A 
comparison  of  this  data  with  that  reported  by  G.  Galvo  [Ref.  14:  p.  61]  showed  that 
there  is  agreement  to  within  the  bounds  of  experimental  uncertainty  except  at  the 
highest  Reynolds  number.  For  comparison,  data  reported  by  Kreith  [Ref.  12:  p.  1253], 
and  Brinich  and  Graham  [Ref.  11:  p.  36j  for  a  convex  surface  were  included  in  the 
graph. 

Using  the  flat  plate  data  as  a  reference  of  comparison,  the  reduction  in  heat 
transfer  on  the  convex  plate  was  manifest  in  a  Nusselt  number  consistently  six  percent 
lower  than  that  of  the  flat  plates  over  the  entire  range  of  How  rates.  This  did  not 
represent  a  large  difference  between  the  convex  and  flat  plate  cases;  however,  six 


percent  was  at  least  one  percent  more  than  the  conservatively  calculated  experimental 
uncertainty  of  the  data.  Therefore,  this  data  implied  that  a  decrease  in  heat  transfer  on 
a  convex  surface  should  be  expected,  though  probably  small. 

A  comparison  of  the  convex  and  concave  surfaces  (see  Figure  4.6)  pointed  out 
that  the  concave  surface  heat  transfer  is  approximately  thirty  percent  greater  than  that 
for  the  convex  surface  at  a  Reynolds  number  of  8000,  and  about  twenty  percent  greater 
at  a  Reynolds  number  of  20,000. 

E.  EFFECT  OF  BUOYANCY 

Tables  3  and  4  include  data  taken  when  the  test  channel  was  oriented  such  that 
free  convection  aided  forced  convection.  Graphical  comparison  for  this  case  is  shown 
in  Figure  4.7. 

Data  taken  for  this  configuration  was  within  the  bounds  of  experimental 
uncertainty,  and  showed  no  pattern  in  differing  from  the  case  where  the  direction  of 
free  convection  opposed  that  of  forced  convection.  These  results  indicated  that  there 
was  no  significant  buoyant  effect  over  the  full  range  of  flow  rates  studied. 
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TABLE  1 

SUMMARY  OF  STRAIGHT  UPPER  WALL  RESULTS 


Rehd 

6ReM  RCM 

8Nuhd  Nu 

6841 

17.30 

±0.0560 

±0.0481 

7856 

19.92 

8588 

21.46 

9327 

23.06 

10168 

24.75 

±0.0576 

±0.0467 

10936 

26.18 

11947 

28.11 

12082 

28.38 

12981 

29.25 

13969 

30.96 

±0.0571 

±0.0485 

14997 

32.66 

15950 

34.62 

17169 

36.85 

17839 

37.30 

±0.0570 

±0.0495 

18927 

39.01 

20008 

40.27 

±0.0570 

±  0.0509 
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TABLE  2 

SUMMARY  OF  STRAIGHT  LOWER  WALL  RESULTS 


Rehd 

Nuhd 

6Re»d  Reh- 

*NuM  Nul 

6071 

15.03 

±0.0581 

±0.0452 

"'065 

17.85 

8116 

19.72 

±0.0573 

±0.0458 

8968 

21.26 

9997 

23.63 

10934 

25.35 

±0.0571 

±0.0464 

12014 

27.17 

12984 

29.11 

14047 

30.89 

±0.0571 

±0.0473 

14989 

32.67 

15981 

34.10 

16984 

35.45 

±0.0570 

±0.0479 

18035 

36.21 

19002 

37.46 

19458 

39.83 

±0.0570 

±0.0474 
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TABLE  3 

SUMMARY  OF  CONCAVE  WALL  RESULTS 


Rehd 

6NuhdNuhd 

5884 

20.04 

±0.0582 

±0.0448 

6025* 

19.98* 

7050 

22.89 

8016 

24.73 

±0.0573 

±0.0450 

8044* 

24.80* 

9094 

26.89 

10025* 

28.08* 

10045 

28.71 

11023 

30.72 

±0.0571 

±0.0451 

11977 

32.72 

13028 

34.90 

13993 

36.19 

±0.0571 

±0.0459 

14982* 

38.05* 

15021 

37.95 

15996 

40.11 

17053 

42.51 

±0.0570 

±0.0464 

17943 

42.90 

19169 

45.26 

20189 

46.80 

±0.0571 

±0.0478 

*  Indicates  measurements  that  were  taken  with  the  test  channel  inverted 
such  that  the  directions  of  free  and  forced  convection  were  coincident. 


TABLE  4 

SUMMARY  OF  CONVEX  WALL  RESULTS 


5  Re.  .  Re.  . 
hd  nd 

*Nuhd'Nuhd 

6065 

15.04 

±0.0581 

±0.0492 

6088* 

15.15* 

7050 

16.85 

S009 

18.87 

±0.0574 

±0.0498 

8034* 

18.88* 

9058 

19.90 

10020 

22.12 

10040* 

22.48* 

11054 

24.32 

±0.0571 

±0.0525 

12003 

24.92 

12988 

27.83 

13944 

29.42 

±0.0571 

±0.0545 

14871* 

30.57* 

15003 

30.72 

15981 

32.20 

17009 

32.84 

±0.0570 

±0.0517 

17898 

35.43 

19015 

36.25 

20210 

37.73 

±0.0571 

±0.0524 

*  Indicates  measurements  that  were  taken  with  the  test  channel  inverted 
such  that  the  directions  of  free  and  forced  convection  were  coincident. 
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Figure  4.2  Straight  Lower  Plate  Results. 
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Figure  4.3  Concave  Plate  Results. 
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Figure  4.6  Comparison  of  Concave.  Convex  and  Straight  Plate  Data 


Figure  4.7  Convex  Concave  Buoyancy  Comparison. 


V.  DISCUSSION 


A.  FLAT  PLATES 

The  experimental  results  for  the  upper  and  lower  straight  sections  of  the  channel 
were  essentially  the  same,  and  are  not  considered  separately  in  this  discussion. 

The  Dittus-Boelter  correlation  for  Nusselt  number  based  on  the  hydraulic 
diameter  of  the  test  channel  and  a  Prandtl  number  of  0.71  yielded  Nusselt  numbers 
thirty  to  thirty-seven  percent  higher  than  those  recorded  in  this  experimentation.  The 
Dittus-Boelter  correlation  is  an  empirical  approximation  that  was  developed  to 
estimate  local  Nusselt  numbers  for  circular  tubes.  The  correlation  itself  is  only  an 
approximation,  and  when  modified  for  non-circular  ducts  where  convection  coefficients 
van-  about  the  periphery,  the  approximation  becomes  less  certain.  Here  the  correlation 
was  used  to  estimate  the  average  Nusselt  number  over  the  length  of  the  heated  section. 
To  a  first  approximation  the  Dittus-Boelter  correlation  can  be  used  to  estimate  the 
average  Nusselt  number  because  the  ratio  of  the  plate  length  to  the  hydraulic  radius 
was  greater  than  ten;  however,  this  assumption  adds  to  the  uncertainty  of  the  resulting 
convection  coefficient. 

The  solution  for  parrallel  plates  of  infinite  extent  developed  by  Kays  and  Leung 
gave  Nusselt  numbers  closer  to  the  experimental  results,  but  was  still  high  by  ten  to 
fifteen  percent,  though  ten  to  fifteen  percent  accuracy  when  estimating  heat  transfer 
coefficients  in  the  turbulent  flow  regime  is  not  unusual.  Even  though  the  straight 
section  of  the  test  channel  approximated  flat  parallel  plates  of  infinite  extent,  some 
sidewall  influence  on  flow  was  expected  and  may  account  for  a  part  of  the  disparity  of 
the  comparison.  The  graphical  comparison  shown  in  Figure  4.1  was  developed  by 
linear  interpolation  of  the  data  from  (Ref.  18:  p.  552),  though  there  was  no  evidence  to 
support  the  assumption  that  the  results  of  the  Kays  and  Leung  solution  arc  a  linear 
function  of  either  Prandtl  or  Reynolds  number;  therefore,  additional  inaccuracy  may 
have  been  introduced  into  the  comparison. 


B.  CURVED  PLATES 


Experimental  and  analytical  investigation  of  flow  and  heat  transfer  on  concave 
and  convex  surfaces  in  a  rectangular  duct  of  high  aspect  ratio  has  thus  far  received 
little  attention;  consequently,  there  is  scant  previous  work  available  for  comparison. 
Some  results  are  available  for  channels  of  modest  aspect  ratios,  but  these  do  not 
generally  compare  well  due  to  the  significant  sidewall  effect  introduced  at  low  aspect 
ratios. 

Studies  of  the  flow  and  heat  transfer  characteristics  in  high  aspect  ratio  channels 
has  been  done  at  the  Naval  Postgraduate  School  (Refs.  9,14,17,19,20,21,22].  Results  of 
this  previous  research  were  comparable  to  the  results  of  this  experimentation  and  were 
highly  depended  upon  to  validate  these  results.  Only  Galyo's  work  is  shown  for 
comparison  in  this  report,  but  his  results  agreed  within  the  limits  of  experimental 
uncertainty  with  those  of  Wilson  [Ref.  17],  and  Daughety  [Ref.  22].  The  major 
improvements  of  this  experimentation  over  that  previously  done  were  the  decreased 
experimental  uncertainty,  especially  at  the  highest  Reynolds  number,  and  the 
preponderance  of  data  that  permitted  empirical  correlations  to  be  determined  to  a  very 
high  level  of  confidence. 

In  1955  Kreith  [Ref.  12:  p.  1253]  investigated  the  influence  of  curvature  on  heat 
transfer  in  a  rectangular  channel  of  aspect  ratio  equal  to  three.  For  a  Prandtl  number 
of  U.71.  graphical  comparison  with  the  results  recorded  in  this  experimentation  showed 
that  Kreiths  results  were  approximately  thirty  percent  lower  (see  Figures  4.3  and  4.5); 
however.  Kreith  s  test  apparatus  not  only  had  a  much  lower  aspect  ratio,  it  was  also 
heated  by  a  heating  dement  that  was  centered  between  the  upper  and  lower  walls  of 
the  channel  instead  of  heating  the  walls  themselves. 

Research  done  in  1977  with  a  channel  of  aspect  ratio  six  by  Brinich  and  Graham 
[Ref.  11:  pp.  35-36]  was  shown  for  comparison  in  Figures  4.3  and  4.5.  The  range  of 
Reynolds  numbers  over  which  they  performed  their  experiments  was  approximately 
$4,000  to  253,000.  For  the  concave  surface.  Brinich  s  and  Graham's  data  extrapolated 
to  a  Reynolds  number  of  15.000  yielded  a  Nusselt  number  about  fifty  percent  higher 
than  that  recorded  in  this  experiment.  Likewise,  an  extrapolation  to  a  RexnolJs 
number  of  15.000  for  the  convex  plate  showed  that  there  was  about  a  thirty  percent 
difference  in  the  results  of  the  two  experiments.  These  differences  are  rather  large,  but. 
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again,  the  two  studies  differed  greatly  in  that  the  aspect  ratios  and  the  range  of 
Reynolds  numbers  investigated  were  much  different. 

Similiar  research  was  done  in  1978  by  L'chida  and  Koizumi  (Ref.  13:  pp. 
1714-1716],  with  channels  of  aspect  ratio  up  to  four;  however,  the  full  report  has  not 
been  translated  into  English.  From  figures,  graphs,  and  equations  included  in  their 
report  it  appears  that  their  test  channel  was  heated  along  its  entire  periphery  and  that 
the  published  results  were  in  the  form  of  an  average  Xusselt  number  for  the 
combination  of  a  convex  and  concave  wall. 

It  is  unlikely  that  an  average  heat  transfer  coefficient  or  Nusselt  number  can  be 
determined  when  both  the  convex  and  concave  walls  are  heated  by  simple 
superposition  of  the  individual  correlations  for  the  convex  and  concave  walls  because 
of  the  complex  nature  of  the  flow  in  a  curved  channel.  Therefore,  there  was  no 
manner  in  which  the  results  detailed  in  Ref.  13  could  justifiably  be  compared  to  the 
results  of  this  experimentation. 

Even  though  the  test  equipment  was  such  that  a  pair  of  heaters  could  be 
energized  to  supply  a  constant  heat  flux  on  opposite  walls  of  the  channel,  no 
conclusive  experiments  could  be  performed  to  determine  an  average  Nusselt  number 
for  this  configuration. 

Referring  to  (eqn.  3.3)  and  denoting  F,  as  the  average  heat  transfer  coefficient  for 
single  plate  heating,  and  Ii2  as  the  average  heat  transfer  coefficient  when  heating  a  pair 
of  opposite  walls,  then: 


«Two  +  Twi>  ^  '  yi 


wh  "  ^blk^ 


(eqn  5. 1 ) 


where  Twh  represents  the  average  temperature  when  one  wall  (inner  or  outer)  is  heated, 
and  the  heat  flux  from  both  the  inner  and  outer  plates  was  assumed  equal  when 
heating  a  pair  of  plates.  The  average  wall  temperature  term  in  the  numerator  of  (eqn. 
5. 1 )  can  be  recast  as: 


(T, 


wo 


Twi>  2 


wh 


oT 


(eqn  5.2  > 
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where  ST  represents  the  difference  between  the  average  temperature  when  heating  a 
pair  of  plates  and  the  single  plate  average  temperature.  Combining  (eqn.  5.1 »  and  (eqn. 
5.2)  yields: 


<T.h  •  Tb  ik  +  «T> 
<Twh  '  Tblk) 

1  +*T(Twh-Th|k) 
1 

1  +  &T  AT 
~  j  — 


(eqn  5.3) 


(Eqn.  5.3)  reveals  that  only  when  the  temperatures  of  both  plates  are  the  same,  such 
that  6T  -  0,  when  heating  a  pair  of  plates,  can  a  true  average  heat  transfer  coefficient 
for  both  walls  be  obtained.  Strictly  speaking  the  temperature  at  all  corresponding 
points  opposite  each  other  on  the  plates  should  be  equal  to  yield  an  average  heat 
transfer  coefficient  for  the  combination  of  both  walls.  The  heaters  used  in  this 
experimentation  were  constant  heat  llux  heaters,  and  temperatures  at  corresponding 
points  on  opposite  walls  could  not  be  matched;  therefore,  conclusive  data  for  the  case 
of  heating  a  pair  of  plates  could  not  be  gathered  using  the  existing  equipment. 

Though  local  heat  transfer  coefficients  and  Nusselt  numbers  were  calculated  at 
the  leading  edge,  at  mid-plate,  and  at  the  trailing  edge  of  the  heated  plates,  no 
conclusive  pattern  was  evident  in  the  data.  The  same  was  true  for  local  temperatures 
More  conclusive  results  could  not  be  obtained  because  there  were  too  few 
thermocouples  spaced  too  far  apart;  that  is.  the  spanwise  thermocouple  spacing  was 
greater  than  the  diameter  of  the  Taylor-Gorrler  vortices  that  formed. 
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VI.  CONCLUSION 


The  results  of  this  experimentation,  detailed  in  Chapter  IV,  indicated  that 
secondary  flow  in  a  curved  rectangular  channel  in  the  form  of  Taylor-Gortler  vortices 
influenced  the  heat  transfer  on  the  concave  and  convex  surfaces  over  a  range  of  flow 
corresponding  to  Reynolds  numbers  from  about  6000  to  20000. 

Correlations  of  Nusselt  number  as  a  function  of  Reynolds  number  for  the  two 
flat  plates,  the  convex  plate,  and  the  concave  plate  were  calculated  to  a  high  level  of 
certainty  and  showed  that  heat  transfer  was  enhanced  by  approximately  twenty  percent 
on  the  concave  surface  as  compared  to  the  Hat  surfaces.  Similiarly,  heat  transfer  over 
the  convex  surface  was  reduced  by  approximately  six  percent  in  comparison  to  the  flat 
surfaces,  f  urthermore,  the  enhancement  or  reduction  in  heat  transfer  on  the  curved 
plates  was  noted  to  be  greatest  at  the  lowest  Reynolds  number  tested,  and  tapered  oil" 
with  increased  How  rate. 

Comparisons  made  of  the  results  obtained  in  this  experimentation  to  that  of 
other  pertinent  research  supported  the  conclusion  that  heat  transfer  is  enhanced  on  a 
concave  surface.  Experiments  performed  lor  the  purpose  of  determining  the  relative 
significance  of  buoyancy  induced  flow  indicated  that  even  at  the  lowest  Reynolds 
number  of  6000,  buoyant  free  convection  was  negligible  compared  to  forced 


VII.  RECOMMENDATIONS 


EQUIPMENT 

1.  Small  oscillations  in  ambient  temperature  (±1°)  noticeably  affected 
temperatures  in  the  channel  and  contributed  to  the  differences  in  the  results  in 
successive  data  taking  runs.  Maintaining  close  control  over  ambient 
temperature  would  make  experimental  results  more  reproducible  and  decrease 
the  time  necessary-  to  gather  data. 

2.  Existing  equipment  allowed  for  testing  only  up  to  a  Reynolds  number  of  about 
20000.  though  the  corresponding  Reynolds  numbers  for  flow  through  a  gas 
turbine  are  much  higher.  A  larger  test  channel  and  higher  capacity  air 
compressor  would  be  required  for  experimentation  at  higher  flow  rates. 

3.  A  great  amount  of  time  was  spent  waiting  approximately  four  hours  for  the  test 
channel  to  reach  steady  state  after  first  energizing  the  heater.  Installing  guard 
heaters  would  significantly  reduce  the  warm  up  time  to  steady  state. 
Additionally,  guard  heaters  would  prevent  heat  loss  through  the  insulation 
surrounding  the  heaters. 

4.  Pressure  sensors  that  could  be  read  directly  by  the  automatic  data  acquisition 
system  several  times  during  the  course  of  a  data  run  would  increase  the 
accuracy  of  flow  rate  measurements  and  permit  successive  data  runs  to  be  taken 
without  the  presence  of  an  operator. 

5.  More  thermocouples  should  be  installed  closer  together  to  improve  the  accuracy 
of  temperature  measurements  and  to  facilitate  studying  the  pattern  of  local 
variations  in  temperature  and  heat  transfer  in  the  spanwise  and  flow-wise 
directions. 

6.  The  single  most  significant  contributor  to  experimental  uncertainty  in  the 
calculated  average  Nusselt  number  and  hydraulic  Reynolds  number  was  the 
uncertainty  in  the  channel  cross  section.  The  uncertainty  in  experimental 
results  could  be  greatly  reduced  if  a  channel  were  constructed  paying  special 
attention  to  the  cross  section  dimensions. 


B.  FURTHER  RESEARCH 

1.  Experimentation  similiar  to  this  but  at  higher  flow  rates  should  be  done  so  that 
a  large  data  base  could  be  established  and  heat  transfer  correlations  could  be 
calculated  for  a  wide  range  of  application. 

2.  Experiments  at  lower  flow  rates  where  buoyancy  is  significant  should  be  done 
so  that  a  suitable  length  parameter  can  be  deduced  for  calculating  the  Grashof 
number  for  curved  channel  flow. 

3.  Different  curvatures  should  be  studied  to  determine  the  effect  of  the  degree  of 
curvature  on  heat  transfer  in  a  curved  channel. 

4.  Thermostatic  control  should  be  included  in  the  experimental  equipment  that 
could  maintain  equal  plate  average  temperatures  when  heating  a  pair  of  heaters 
on  opposite  walls  so  that  experimental  results  for  channel  average  heat  transfer 
coefficients  could  be  gathered. 
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APPENDIX  A 

DATA  REDUCTION  PROGRAM 


1000  1  FILE  NAME:  MAIN 

1010  !  LAST  REVISION  DATE:  19  OCT  1986 

1020  !  EDITED  BY:  LT.  J.R.  HAWK  III 

1030! 

1040  !  PROGRAM  FOR  GATHERING  AND  REDUCING  DATA  WHEN  HEATING  ONE  COPPER  PLATE, 

10S0  !  WITH  S  MINUTES  BETWEEN  RUNS.  PROGRAM  IS  MODIFIED  TO  ACCOUNT  FOR 

1060  I  BROKEN  THERMOCOUPLE  CHANNELS  S,  t6,  2t,  AND  26  WHICH  EFFECT  OUTER 
1070  !  CURVED  HEATER.  THERMOCOUPLE  CHANNEL  47  WHICH  EFFECTS  INNER  CURVED 

1080  !  HEATER,  AND  THERMOCOUPLE  CHANNEL  11  WHICH  EFFECTS  OUTER  STRAIGHT  HEATER 

1090! 

1100  COM  /Co/  D(7) ,Aa( 76,2) 

1110  DIM  Emf ( 79) , T( 79) 

1120! 

1130  1  CORRELATION  FACTORS  FOR  EMF  TO  DEGREES  CELSIUS  FOR  CU-CO  T/C. 

1140! 

1 1  SO  DATA  0. 10086091 ,25727.94369.-767345.8295.78025595.81 
1160  DATA  9247486589, 6. 97688E*11. -2. 66192EM3, 3. 9407BEH4 
1170  READ  DC) 

1180  PRINTER  IS  701 
1190  BEEP 

1200  INPUT  "ENTER  MONTH,  DAY,  AND  TIME  (MM*.DD:HH:MH:SS> *  ,  T  ime* 

1210  OUTPUT  709; " TD" ; T lme$ 

1220  BEEP 

1230  INPUT  "ENTER  INPUT  MOOE  (0-3054A,  1-FILE)".  Ii* 

1240  IF  Im-1  THEN 
1250  BEEP 

1260  INPU1  "ENTER  NAME  OF  EXISTING  DATA  FILE" ,01df i le* 

1270  PRINT  USING  "20X, ""These  Results  Are  From  Data  File:  * " , 1 OA" ;01df i leS 

1280  PRINT 

1290  PRINT 

1300  ELSE 

1310  BEEP 

1320  INPUT  "ENTER  NAME  OF  NEW  DATA  FILE?" ,Neuf 1 le$ 

1330  PRINT  USING  *20X,""Data  For  This  Run  Is  Stored  In  Data  File:  "“.10ft  ;Neuf 
l  le<S 

1340  PRINT 
1350  PRINT 
1360  ENT)  IF 
1370! 

1380  !  ENTER  THE  CHANNEL  CONFIGURATION. 

1400*  PRINT  USING  "9X,*"THE  FOLLOWING  CHANNEL  CONFIGURATION  WAS  SELECTED:""" 

1410  PRINT 

1420  BEEP  „  . 

1430  INPUT  “SELECT  HEATER  TYPE  TO-STRAIGHT,  1 -CURVED  1 ",  I  type 

1440  IF  Itype-0  THEN  PRINT  USING  "12X, ""Heating  STRAIGHT  Test  Section.1 

1450  IF  I  type- 1  THEN  PRINT  USING  " 12X, * "Hea t i ng  CURVED  Test  Section. """ 

1460  BEEP  .  , 

1470  INPUT  "SELECT  HEATER  POSITION  (O-OUTER,  1 -INNER) ", Ipos 
1480  If  Ipos-0  THEN  PRINT  USING  " 12X ,* "Hea t t ng  OUTER  Plate."" 

1490  IF  Ipos-1  THEN  PRINT  USING  "12X, ""Heating  INNER  Plate.""" 

1500  PRINT 
1510! 
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1520 
1530! 
15*40 
1550 
1560 
1570 
1580 
1590 
1600! 
1610 
1620! 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740! 
1750 
1760! 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1360 
1870 
1880 
1890! 
1900 
1910! 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
dd  i 
2000 
dd) 

20  1  0 
s  ) " "  “ 
2020 
s  )  “  “ " 
2030 
2040 
2050 
2060! 
2070 
2080! 


!  ENTERS  FILES  WITH  T/C  CALIBRATION  COEFFICIENTS  BASED  ON  HEATER  TYPE. 

IF  Itype-0  THEN  ASSIGN  «Filec  TO  "STRA_CQEF* 

IF  Itype-1  THEN  ASSIGN  SFiiec  TO  *CURV_COEF“ 

FOR  1-0  TO  76 

ENTER  $F  i  led Aa(  1 , 0 ) , Aa(  1 , 1  > . Aa(  1 ,2) 

NEXT  I 

ASSIGN  dFilec  TO  * 

!  INSERT  DATA  STORAGE  DISK. 

BEEP 

INPUT  "CHANGE  DISK  AND  HIT  ENTER". Ok 

IF  Im-1  THEN  ASSIGN  SFile  TO  OldfileS 

IF  Im-0  THEN 

CREATE  BDAT  NewfileS.40 

ASSIGN  9F  i  le  TO  NeufileS 

END  IF 

ENTER  709: Time* 

CLEAR  709 

PRINT  USING  "20X, ""Month.  Day.  and  Time:  — , 1 4A" ; Time* 

PRINT 

!  FOR  DATA  RUN  READ  AND  STORE  ALL  RAW  EMF  VALUES. 


IF  Im-0  THEN 

OUTPUT  709: "AR  AFOO  AL79" 
OUTPUT  722 ; *F  1  R1  T1  Z1  FL1 " 
FOR  1-0  TO  79 
OUTPUT  709; "AS" 

ENTER  722;Emf ( I) 

NEXT  I 

OUTPUT  ®F i le;Emf (*) 

ELSE 

ENTER  8Ftle:Emf<*) 

END  IF 

OUTPUT  709; "TD" 


!  CONVERT  AND  CALIBRATE  ALL  T/C  READINGS  TO  DEGREES  CELSIUS. 

PRINT  USING  "9X, ""THE  FOLLOWING  DATA  WERE  RECORDED:""" 

PRINT 

FOR  1-0  TO  79 

IF  1-38  OR  1-39  OR  1-77  OR  1-78  THEN 

IF  ( Ipos-0  AND  1-38)  OR  (Ipos-1  AND  1-77)  THEN 

Iadd-0 


IF  Ipos-1  THEN  Iadd-39 


1 tage  ( Vh ) 

-  "“.2D. 

3D . "  “ 

(V) 

“ * " :Emf ( 38 ♦  la 

Voltage  i V r ) 

-  "".2D. 

3D ,  “  “ 

l  V  > 

“““ : Emf ( 39*  la 

12X,  '"'Precision 

Resistor 

( Rpr ) 

- 

2.01729  (Ohm 

12X. "“Precision 

Resistor 

( Rpr ) 

* 

2.00839  (Ohm 

PRINT 
END  IF 
ELSE 


!  IF  T/C  OUT  OF  RANGE  PRINT  A  OUE . 
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2090  Aem*  -Emf ( I ) 

2100  IF  Aemf < 7 . 09E-4  OR  Aemf >3.630E-3  THEN 
2110  T( I ) *999.99 
2120  ELSE 

2130  CALL  Tvsv(Emf ( I ) . Tt ) 

21  AO  IF  1-79  THEN 

2150  T< I ) - Tt 

2160  GOTO  2210 

2170  END  IF 

2180  T ( I ) -FNTcorr ( T t , I ) 

2190  END  IF 

2200  END  IF 

2210  NEXT  I 

2220! 

2230  !  PRINT  ALL  T/C  READINGS. 

2240! 

2250  PRINT  USING  "12X, "OUTER  PLATE:"* 

2260  PRINT  USING  *14X,"T/C  Number:  1  2  3  4  S 

gH  It  • 

2270  PRINT  USING  *14X,"Temp  (C):  *  *  .  6(  3D.  DD.2X)  *  J  T(  0 ) ,  T(  1 ) .  T(2) .  T<  3  > .  T<  4) .  T 

(5) 

2280  PRINT  USING  *14X."T/C  Number:  7  8  9  10  11 

12"" 

2290  PRINT  USING  *14X,"Temp  (Cl:  "  .6<3D.DD.2X)*;T(6> ,T(7) ,T(8) ,T(9) ,T( 10) , 

T<  1 1 ) 

2300  PRINT  USING  *14X,"T/C  Number:  13  14  15  16  17 

1 8“ " 

2310  PRINT  USING  *14X,"Temp  (C):  *  *  ,6(  3D.DD  ,2X )  *  :  T  ( 12) .  T(  1 3  > .  T(  1 4 ) ,  T<  15  > .  T( 

16) ,T<17) 

2320  PRINT  USING  *14X,“*T/C  Number:  t9  20  21  22  23 

24"“ 

2330  PRINT  USING  "14X."Temp  <C):  ”  ,6(  3D.  DD.2X )  *  ;  T(  1 8) .  T<  1 9) .  T(20 ) ,  T  ( 21  > ,  T( 

22 )  T ( 23 ) 

2340  PRINT  USING  *14X,"T/C  Number:  25  26  27  28  29 

3C"“ 

2350  PRINT  USING  *14X,"Temp  (C) :  "  ,6(  3D. DD.2X )  *  ;  T(24 ) .  T< 25 ) .  T ( 26 ) ,  T< 27 ) .  T ( 

28) . Tt 29) 

2360  PRINT 
2370! 

2380  !  CALCULATE  AVERAGE  TEMPERATURE  OF  OUTER  HEATED  WALL. 

2390! 

2400  Sum-0. 

2410  IF  I  type-0  THEN 

2420  FOR  1-0  TO  29 

2421  IF  1-11  THEN  2440 
2430  Sum-Sum+T ( I ) 

2440  NEXT  I 

2450  Tuo-Sum/29 
2460 ! 

2470  !  CORRECTION  FOR  BROKEN  T/C  CHANNELS  5.  16.  21.  AND  26  IN  OUTER  CURVED  HE 

ATER. 

2480! 

2490  ELSE 

2500  FOR  1*0  TO  29 

2510  IF  1-5  OR  1-16  OR  1-21  OR  1-26  THEN  2530 
2520  Sum-Sum+  T ( I ) 

2530  NEXT  I 
2540  Two-Sum/26 
2550  END  IF 

2560  PRINT  USING  "12X. "OUTER  INSULATION:"" 


2570  PRINT  USING  *14X,"*T/C  Number 


31  t  T insol ) 


32  *Tmso2)  33  (Tinso3> 


j'l’j 


M  « 

2580  PRINT  USING  "14X,"‘Temp  1C):  "" .  3(  3D.0D.8X )  * ;  TOO) ,  T(  31 ) ,  T<  32) 

2590  T  i  nsol *  T ( 30  > 

2600  T  mso2*T (31  ) 

2610  T  inso3*T (32) 

2620  PRINT 

2630  PRINT  USING  * 12X, ""ORIFICE  TEMP  (Torf)  -  "" .3D.DD, "*  (C)*"";T(33) 

2640  PRINT 

2650  PRINT  USING  *12X.**INLET  TEMPERATURE :" "* 

2660  PRINT  USING  "14X,*"T/C  Number:  35  36  37  38“"* 

2670  PRINT  USING  "14X,""Temp  (C):  " " , 4< 3D. DD . 2X ) " ; T( 34 ) , T ( 35 ) . T( 36 ) , T< 37 ) 

2680  PRINT 

2690! 

2700  !  CALCULATE  AVERAGE  INLET  TEMPERATURE. 

2710! 

2720  Sum-0. 

2730  FOR  1-34  TO  37 
2740  Sum-Sum-T(I) 

2750  NEXT  I 
2760  T  in-Sum/4 

2770  PRINT  USING  "12X, ""INNER  PLATE:""" 

2780  PRINT  USING  "14X.*"T/C  Number:  41  42  43  44  45 

46""* 

2790  PRINT  USING  "14X.*"Temp  <C):  "" ,6< 3D.DD.2X) * ; T< 40 ) . T ( 4 1 ) . T< 42 ) . T l 43) , T ( 

44) ,T(45> 

2800  PRINT  USING  "14X."“T/C  Number:  47  48  49  50  51 

£  -1  •»  «  « 

2810  PRINT  USING  "14X,*"Temp  <C>*.  *"  .613D.DD.2X) * : T ( 46) . T (  47 ) . T ( 48 > .  T ( 49 > .  T< 

50)  ,T(51  ) 

2820  PRINT  USING  "14X,""T/C  Number:  53  54  55  56  57 

53-.. 

2830  PRINT  USING  "14X.""Temp  (C>:  " * ,6( 3D.DD .2X ) " ; T ( 52 ) . T(53 > . T ( 54 ) , T ( 55 ) . T ( 

56) . T< 57) 

2840  PRINT  USING  "14X,““T/C  Number:  59  60  61  62  63 

64“““ 

2850  PRINT  USING  *14X.""Temp  <C>:  " " ,6( 3D. DD.2X ) " : T( 58 ) . T ( 59 ) . T ( 60 ) . T ( 61 ) . T< 

62)  . T ( 63 ) 

2860  PRINT  USING  "14X,""T/C  Number:  65  66  67  68  69 

70“““ 

2870  PRINT  USING  *14X,""Temp  <C):  " " ,6( 3D. DD.2X ) " ; T< 64 ) . T 1 65  ) , T < 66  ) , T ( 67  )  , T ( 

63 )  . T ( 69 ) 

2880  PRINT 
2890! 

2900  !  CALCULATE  AVERAGE  TEMPERATURE  OF  INNER  HEATED  NALL. 

2910! 

2 920  Sum-0. 

2930  IF  I  type-0  THEN 
29ui)  FOR  I*mi)  TO  69 
2950  Sun-Sum-TlI) 

2960  NEXT  I 
2970  Tiui-Sum/30 
2980! 

2990  !  CORRECTION  FOR  BROKEN  T/C  CHANNEL  47  IN  INNER  CURVED  HEATER. 

3000! 

3010  ELSE 

3020  FOR  1-40  TO  69 
3030  IF  1-47  THEN  3050 
3040  Sum-Sum- T (  I ) 

3050  NEXT  I 


V  V 
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3060  Twi-Sum/29 
3070  END  IF 

3080  PRINT  USING  * 12X. **  INNER  INSULATION: “** 

3090  PRINT  USING  *  1 4X T/C  Number :  71  ( T i ns 1 1 )  72  (Tinst2)  73  (Tinsi3> 

•*  • 


3100  PRINT  USING  "14X,*"Temp  <C):  " " ,3(30.DD, 8X) * : T( 70 ) , T< 71 1 . T( 72) 

3110  T  i  ns 1 1 *T ( 70 ) 

3120  T insi2“T<  7 1 ) 

3130  T  insi3*T ( 72) 

3140  PRINT 

3 ISO  PRINT  USING  *  1 2X ." "OUTLET  TEMPERATURE:*” 

3160  PRINT  USING  *14X,”T/C  Number:  74  75  76  77”* 

3170  PRINT  USING  "14X.  “Temp  (C):  ”  , 4( 3D . DD.2X ) * ; T( 73 ) , T( 74 ) , T ( 75 ) , T ( 76) 

3180  PRINT  USING 

3190! 

3200  !  CALCULATE  AVERAGE  OUTLET  TEMPERATURE. 

3210! 

3220  Sum-0. 

3230  FOR  1-73  TO  76 
3240  Sum-Sum*T ( I ) 

3250  NEXT  I 
3260  Tout-Sum/4 
3270! 

3280  !  CALCULATE  DIFFERENCE  BETWEEN  CHANNEL  INLET  AND  OUTLET  TEMPERATURES. 

3290! 

3300  Tdif f-Tout-Tin 
3310! 

3320  !  CALCULATE  CHANNEL  BULK  TEMPERATURE. 

3330! 


3340  Tblk*( T i n+Tout ) /2 

3350! 

3360  !  PRINT  CALCULATED  TEMPERATURES. 

3370! 

3380  PRINT  USING  *9X,”THE  FOLLOWING  TEMPERATURES  WERE  CALCULATED:”* 


3390  PRINT 

3400  PRINT  USING  *  1 2X 'Average  Outer  Wall  Temperature  (Two)  -  ”,3D 

DD.”  (  C )  "  *  *  ;  Two 

3410  PRINT  USING  “ 1 2X, " “Average  Inner  Wall  Temperature  (Twi)  -  ”,3D 

DD.”  (C)”";Tui 

3420  PRINT  USING  “ 1 2X , " “Averaae  Outlet  Temperature  (Tout)  *  ”,3D. 

DD.*’  ( C ) - : Too t 

3430  PRINT  USING  " 1 2X , "“Average  Inlet  Temperature  (Tin)  -  "*,3D. 

DD.”  (  C  ) - :  T  i  n 

3440  PRINT  USING  * 12X, * "Channe  l  Inlet  and  Outlet  Temp  Diff  (Tdiff)  ■  ”.3D. 

DD.”  Tdiff 

3450  PRINT  USING  " 12X.  ""Average  Bulk  Temperature  (Tblk)  -  ”,3D. 

nn  -  *  i i  *>.TkiL 


3460  1 

34  7,)  !  CALCULATE  DELTA  T  (Tdel!  FOR  USE  IN  CALCULATING  H  LATER. 

34'iO  1 


,4.111 

Tde lo- Two- Tb  lk 

3500 

Tdel i*Twi-Tblk 

35  1  0 

IF  Ipos-0  THEN  Tdel-Tdelo 

3520 

IF  Ipos-1  THEN  Tdel-Tdeli 

3550 

PRINT  USING  "12X."“Mean  Temperature  Di 

f  f  erence 

(Tdel) 

DD.” 

( C  )"""; Tde 1 

3540! 

3550 

!  CALCULATE  DIFFERENCE  BETWEEN  OUTER 

AND  INNER 

WALL. 

3560' 

5570 

T wd  1 1  ('  •  ABS  (  Two  -  Tw  i  ) 

53 


3580  PRINT  USING  *12X,*’0uter  and  Inner  Mall  Temp  Difference  iTudiff)  *  ' 

DD.""  (C)""";Twdiff 

3590! 

3600  !  CALCULATE  LOCAL  CHANNEL  TEMPERATURES  AT  INLET.  MIDDLE  AND  OUTLET. 

3610! 

3620  Tcin-T in*( ,5*Tdi f f / 1 2 ) 

3630  Tcaid* T i n*< 4 .5* Td i f f /1 2 > 

3640  lcout*Tin*(  1 1 . 5«Tdif  f  /1 2) 

3650  PRINT  USING  *  12X ," "Average  Local  Channel  Temp  at  Inlet  (Tcin)  *  ' 
DD.""  (C) — ; Tcin 

3660  PRINT  USING  * I2X. ""Average  Local  Channel  Temp  at  Middle  (Tcmid)  •  ' 
DD .  “ "  (C)*"":Tc«ud 

3670  PRINT  USING  “ 12X ,* "Average  Local  Channel  Temp  at  Outlet  (Tcout)  *  1 

DD , ‘ "  (C)“"‘:Tcout 

3680! 

3690  !  CALCULATE  AVERAGE  LOCAL  HEATED  HALL  TEMPERATURE  AT  INLET. 

3700! 

3710  Sum-0. 

3720  Iadd-0 

3730  IF  Ipos-I  THEN  Iadd-40 
3740! 

3750  !  CORRECTION  FOR  BROKEN  T/C  CHANNELS  5.  AND  47  CURVED  HEATER  ONLY. 

3760! 

3770  IF  I  type- 1  THEN 
3780  FOR  1-0  TO  6 

3790  IF  (Ipos-0  AND  1-5)  OR  (Ipos-1  AND  1-6)  THEN  38)0 
3800  Sum-Sum*T < I ♦ Iadd ) 

3810  NEXT  I 
3820  Tum-Sun/6 
3830  ELSE 
3840  FOR  1-0  TO  7 

3850  IF  l Ipos-0  AND  1-7)  OR  (Ipos-I  AND  1-6)  THEN  3870 
3860  Sum* Sum* T ( I  *Iadd  > 

3870  NEXT  I 
3880  Tum-Sum/7 
3890  END  IF 

3900  PRINT  USING  *12X.""Avg  Local  Heated  Hall  Temp  at  Inlet  (T«m)  -  1 

DD.""  (C) " " * ; Twin 

3910! 

3920  !  CALCULATE  AVERAGE  LOCAL  HEATED  HALL  TEMPERATURE  NEAR  MIDDLE. 

3930! 

3940  Sum-0. 


3960  !  CORRECTION  FOR  BROKEN  T/C  CHANNEL  T6  CURVED  HEATER  CNl T .  AND  T/C 

3961 (CHANNEL  11  STRAIGHT  HEATER  ONLY. 

3970! 

3980  IF  I  type-!  THEN 
3990  FOR  1-9  TO  16 

4000  IF  (Ipos-0  AND  1*10)  OR  1 Ipos-0  ANT)  r-!6i  OR  <Ipos*1  AND  I-Ai  'hE’i  -iO.  . 
4010  Sum-Sum*T( I  * Iadd > 

4020  NEXT  I 

4030  IF  Ipos-0  THEN  Turn id-Sum/6 . 0 

4031  IF  Ipos-1  THEN  Twmid*3um/7. 3 
4040  ELSE 

4050  FOR  1-9  TO  16 

4060  IF  (Ipos-0  AND  1-10)  OR  (Ipos*i  AND  1-9)  OR  (Ipos-0  AND  !-<li  ’MfN  aoeo 
4070  Sum-Sun* T ( I  * Iadd ) 

4080  NEXT  I 

4081  IF  Ipos-0  THEN  Tumid-Sum/6.0 
4090  IF  Ipos-1  THEN  Turn id-Sum/7 . 0 


4100  END  IF 

4110  PRINT  USING  “  1 2X , “  "Avg  Local  Heated  Wall  Temp  Near  Middle  (Tumid)  •  “,3D 

DO.”  «C) - ;  Tumid 

4120! 

4130  !  CALCULATE  AVERAGE  LOCAL  HEATED  WALL  TEMPERATURE  AT  OUTLET. 

4140! 

4150  Sum-0. 

4160! 

4170  !  CORRECTION  FOR  BROKEN  T/C  CHANNEL  26  CURVED  HEATER  ONLY. 

4180! 

4190  IF  I  type- 1  AND  Ipos-0  THEN 
4200  FOR  1-23  TO  29 
4210  IF  1-26  THEN  4230 
4220  Sum-Sum+T ( I ) 

4230  NEXT  I 
4240  Tuout-Sum/6 
4250  ELSE 

4260  FOR  1-23  TO  29 
4270  Sum-Sum+T( I+Iadd) 

4280  NEXT  I 
4290  Tuout-Sum/7 
4300  END  IF 

4310  PRINT  USING  " 12X, "*Avg  Local  Heated  Wall  Temp  at  Outlet  (Tuout)  •  “‘.3D 

DD . '  *  (C>  "“I  Twout 

4320  PRINT 

4330  PRINT 

4340! 

43S0  !  ENTER  ORIFICE  PRESSURE  DATA. 

4360! 

4370  BEEP 

4380  INPUT  ’ENTER  PATH  (inHg).  DPM(inH20>.  P1M(inH20).  RE1* .Patm.Dpm. Pirn. Ret 
4390  PRINT  USING  -9X.--THE  FOLLOWING  ORIFICE  DATA  WERE  ENTERED: ■“* 

4400  PRINT 

4410  PRINT  USING  -12X.--Patm(  inHg)  DPM(inH20)  P1M(inH20>  REr" 

4420  PRINT  USING  *  1 4X ,2(2D. 00 , 1 0X ) , (2D. DD . 7X ) ,5D . D* ;Patm, Dpm.P1m.Re1 
4430  PRINT 
4440! 

4450  !  CONVERT  PRESSURE  READINGS  TO  SI  UNITS  AND  P1M  TO  ABSOLUTE. 

4460! 

4470  Patm-Patm*3376.8 
4480  Pde 1 -0pm*248 . 84 
4490  PI -Pa tm- (P1m*248.84 ) 

4600  Tor»-T ( 

4510! 

4620  !  PHYSICAL  PROPERTIES  AND  CONSTANTS. 

4630' 

4640  R-286.987  !F0R  AIR. 

4550  Rho-P1/<R*( Tort »273. ISM  !  DENSITY  OF  AIR  AT  ORIF  PLATE. 

4C.  nu  Gamma- 1.40  'FOR  AIR, 

46 7 0  Co-1006  'Cp  GOOD  FOR  Tblk  BETWEEN  12  AND  33  (DEG  C). 

46.JIJ  be  -I.u 

46  }<l  Sigma-5. 663E-8 

46.10  Ecu-. 12  1EMISSIVITT  OF  COPPER. 

4tl  1  0  1 

4620  !  DIMENSIONS  OF  CHANNEL. 

4610' 

46uO  Dc-. 006350 
4650  Hid*. 2640 
4660  R  i  • . 297 
4t-,  7 0  Puet  -2»  •  Dc -W  id  ) 

4ho0  Ac-Qc-Wid 
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4690  Dhd*4*Ac/Puet 
4700  Apl-. 0774192 
4710! 

4720  !  PROPERTIES  OF  INSULATION. 

4730! 

4740  Xms-.0127 

4750  Kins-3.8E-2 

4760! 

4770  !  ENTER  ORIFICE  CONFIGURATION. 

4780! 

4790  BEEP 

4800  INPUT  ‘SELECT  DIAMETER  OF  ORIFICE  10-.S334,  1-1.0755  ( inches) ) *  ,Si  10 
4810  IF  S 1 ze-0  THEN  Dorf-. 013548 
4820  IF  Size-1  THEN  Dorf-. 027318 
4830! 

4840  !  DIMENSIONS  OF  PIPE  AND  ORIFICE  PLATE. 

4850! 

4860  A-(PI*Dorf ‘2)/4 
4870  Dpipe-. 051772 
4880  Apipe-PI»Dpip«‘2/4 
4890  Beta-Dor f /Dp  ip* 

4900  PRINT  USING  *  1 2X.  "Dorf  <  •  >  A f » ’ 2 )  Beta*” 

4910  PRINT  USING  *  1 2X , ( Z. 6D ,6X ) , ( Z . 3DE.6X ) , Z . 4D* ; Dorf , A , Be ta 
4920  PRINT 
4930  PRINT 
4940! 

4950  !  CORRELATION  FOR  EXPANSION  FACTOR  BASED  ON  ID  AND  1/2D  TAPS. 

4960! 

4970  Y*  1  -  ( .4l  +  .3S*Beta‘4)*(Pdel/( Gae»a*P1 > ) 

4360! 

4990  !  CALCULATION  OF  FLOW  COEFFICIENT  (K)  IN  SI  UNITS. 

5000! 

5010  B-. 0002*2. 794E*5/0pip«M .0038* 1.0 16E-5/Dpipe>*< Beta *2M 16. 5*1. 968504E*2*Dp 

ipe)»Beta‘ 16) 

5020! 

5030  !  K1  AND  K2  ARE  DUMMY  VARIABLES  TO  CALCULATE  Ko. 

5040! 

5050  K1-. 6014-5. 3974064E-3/Dpipe-. 25 

5060  K2-( .376»2.8971138E-2/Dpipe*.25)*M .6129E-7/(Dp»p#‘2*B*ta‘2*6.35E-5«0pipe) 

♦Beta "4+1 . 5 • Beta' 16) 

5070  Ko-K 1 *K2 

5080  K-K0*1 000*B/R#1 * .5 

5090! 

5100  !  Mu  CORRELATION  GOOD  FOR  Torf  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

5110! 

5120  Mu-4.6971429E-8*Torf ♦ 1 . 7 I94722E-5 

51  30! 

5140  !  CALCULATE  Mdot  AND  REpip*  AND  COMPARE  TO  PREDICTED  RE  I . 

5 1  5'j  1 


5160  Ado  t -K»A»  Y ♦ i 2-Oc  *Rho*Pd# 1 ! ‘ .5 
5170  Rep  ip# *  1 Mdo  t  »0p ip* ) / < Mu-Ap  tp# 1 
5180  D 1 ) f • <  A8S<  Re  I  -R#p  ip* ) /R# I >♦100 
5130  Rel -Repip* 

5200  IF  0  i  f  f  > . 00 1  THEN  5080 

5210  PRINT  USING  *9X,"THE  FOLLOWING  DATA  WERE  CALCULATED: 
5220  PRINT 

5230  PRINT  USING  * 12X.  "Or  1 f ice  Expansion  Factor  <Y> 

5240  PRINT  USING  ‘  1 2X  .  *  "Or  1  f ice  Plow  Coefficient  (M 
5250  PRINT  USING  *  1 2X .  "D*nS 1 1 y  Based  on  T0rf  (Rho) 
ji"“  :Rho 


. Z . 4D. * :  y 
,Z.4D.“:k 

, Z . 4D .  "  <  k  a/w  4 
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S26Q  PRINT  USING  * 1 2X, “ 'Viscosi ty  Based  on  Torf  (Hu)  ■  ".2Z.3DE,"  <kg/m.s> 
5270 ! 

5280  !  Hu  CORRELATION  G000  FOR  Tblk  BETWEEN  17  AND  41  DEGREES  CELSTUS. 

5290! 

5300  Mu*4.697f 429£-8*Tblk  *1 . 7 194722E-5 
5310  Red* ( Hdot*Dc I / ( Hu*Ac ) 

5320  Rehd* < Hdot »Dhd) / ( Hu*Ac) 

5330! 

5340  !  CALCULATION  OF  POWER  INTO  HEATER  PLATE. 

5350! 

5360  Iadd-0 

5370  IF  Ipos-1  THEN  Iadd-39 

5380  IF  Ipos-0  THEN  Rpr-2. 01729 

5390  IF  tpos-1  THEN  Rpr-2. 00839 

5400  Qp*( E»f ( 38*Iadd) »Enf ( 39+ I add) ) /Rpr 

5410! 

5420  !  CALCULATION  OF  HEAT  INPUT  TO  AIR. 

5430! 

5440  Qat r*Hdot*Cp*< Tout-Tin) 

5450! 

5460  !  CALCULATE  LOCAL  HEAT  TRANSFER  COEF .  AND  NUSSELT  NUHBER . 

5470! 

5480  ((air- 7. 7257 l43£-S*Tcin*. 024 165836 
5490  Hm-Qair/(Apl*(  Tujn-Tcin) ) 

5500  Nuin-(Hin*Dhd)/Kair 

5510  Ka  i  r-7 . 7257 1 43E-S*Tc«id* . 024165836 

5520  H«id*Qair/(Apl*( Turn  id- Tew  id) > 

5530  Nu*ud-<H«ud*Dhd)/Kair 

5540  Kair-7. 72571 43E -5* Tcout* .024 165836 

5550  Hout-Qa i r /( Ap  l»( Tuout- Tcou  t ) ) 

5560  Nuout-(Hout»Dhd)/Kair 
5570! 

5580  !  Katr  CORRELATION  GOOD  FOR  Tblk  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

5590! 

5600  Ka  ir- 7. 7257 143E-5*Tb Ik *.0241 65836 
5610! 

5620  !  CALCULATE  AVERAGE  HEAT  TRANSFER  COEF.  AND  NUSSELT  NUHBER. 

5630! 

5640  Havg*Qair/(Apl»Tdel ) 

5650  Nuavg-(Havg»Dhd)/Kair 
5660 ! 

5670  !  CALCULATE  HEAT  LOSSES. 

5630 ! 

56  30  Qlol  •  ( K  i r»s*Ap  1  •  (  T  inso  I  - T  i  nso2  )  1  7X i ns 

5700  Q lo2* i K ins*Ap 1 • l T  mso2- T inso3  > I /X i ns 
5710  0  lo- <  Q lo 1 +Q lo2 ) /2 

5720  01  1 1  •  t  K  i  r>s»Apl  •  (  T  ms  1 1  -  T  ms  i2 )  i/Kms 

57  30  0 1  i2*  <  K  i ns*Ap  1 »(  T  msi2  - T i ns 1 3 ) )  /  '( i  ns 

5740  01  i-(Ql il *01 12)/2 

5750  R  r  *2 . 0  *  ( (  I-£cu)/<Apl-*£cu)  )  *  <  1  . 0/Hp  l  > 

5760  Qr  *  <  S i qma* ( A8S( ( Two*273,1S)‘4-( Tui*273.lS) '41 i )/Rr 

5770  Ode i *0p-0a i r-Qlo! -Or -01 1 1 

5780! 

5790  !  CALCULATE  DEAN  NUMBER. 

SHOO! 

5810  De* .S»Red* ( ( Dc/2 ) /R 1 1  *  .5 

5320  PRINT  USING  "  1 2X .  *  *V  iscos  i  ty  Based  on  Tblk  (  Hu  )  -  ".2Z.3DE."  <kq/m.s) 

4  *  •  i  y 

,  flu 

5330  PRINT  USING  "ICX, "Therm  Cond  0ased  on  Tblk  <Kair>  -  ".Z.3DE."  <W/m.*l* 
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*■ ;Kair 
5840  PRINT 

5850  PRINT  USING  *12X.“*Patn 
5860  PRINT  USING  *12X.**Pdal 
5870  PRINT  USING  *12X.**P1 
5880  PRINT  USING  *12X,**Mdot 
5890  PRINT  USING  *  1 2X . * *Rap ipa 
5900  PRINT  USING  *12X.**Rad 
5910  PRINT  USING  *12X.**Rehd 
5920! 

5930 
5940! 

5950 
5960 
5970 

atts)**’  ‘.Qp.Qair 
5980  PRIN1  USING  *12X.**Qlol 
(Ha t  ts) * “ “ :Qlo1 .01 1 1 
5990  PRINT  USING  *12X,**Qlo2 
(Watts)”";01o2,Qli2 
6000  PRINT  USING  *12X.**Qlo 
(Ha t ts i “  *  * iQlo.Ql t 
6010  PRINT  USING  *12X,**Qr 
(Hatts)"*“:Qr.Qdel 
6020  PRINT 

6030  PRINT  USING  *12X.**Hin 
in.Nuin 

6040  PRINT  USING  *12X,**Hmd 
mid. Nun  id 

6050  PRINT  USING  *12X.**Hout 
out .Nuout 


** .6D.DD. **  (N/n‘2) *** ;Patm 
*■,40.00,**  ( N/n ' 2  >  *  *  * ;Pd« i 
**,6D.DD. **  (N/»-2)***:Pl 
*  “ ,  Z .  4D .  *  *  (kg/s)**'  Jfldot 
** ,6D.D* : Rep  ip* 

** ,60.0* :R«d 
**,6D.D“;R«hd 


30. 

3D.** 

(Hatts)**. 

7X ,  * 

*Qair  - 

.2D 

.30.** 

(Hatts)** 

,7X. 

*  *  0 1 1 1 

.20 

.30.** 

(Hatts)** 

.  7X, 

**01i2 

.20 

.3D.** 

(Hatts)** 

,7X. 

**Q1 1 

.20 

.30.** 

(Hatts)** 

,7X. 

**Qdsl 

3D. 

DO.** 

(U/*'2C) ** 

.  7X , 

**Nuin 

3D. 

DD.** 

(H/n'2C) ** 

.  7X , 

**Nun.d 

3D. 

DD.** 

(U/n'2C) ** 

,  7X , 

**Nuout 

!  ONLY  PRINT  DEAN  NUMBER  UHEN  HEATING  CURVED  SECTION. 


IF  Itype-1  THEN  PRINT  USING  *12X.**0a 
PRINT 

PRINT  USING  “12X. “*Qp 


**.40.D*;0a 


6060  PRINT  USING  *12X,**Havg  -  •*,30.DD,**  (H/»*2C) ** , 7X , * *Nuavg  • 

avg.Nuavg 

6070  PRINT 

6080  PRINT 

6090  BEEP 

6100  INPUT  "HILL  THERE  BE  ANOTHER  RUN?  ( 1 • Y, 0-N) * ,Go_on 
6110  IF  Go.onOO  THEN 
6120  PRINT  USING  *9.#* 

6130  IF  In-0  THEN  HAIT  300 
6140  GOTO  1700 
6150  END  IF 

6160  PRINT  USING  *8X,**EN0  OF  RUN*** 

6170  PRINT  USING  *$.»* 

6180  PRINTER  IS  1 
6190  ASSIGN  sFile  TO  * 

6200  END 


6210! 

6220  !  CONVERTS  EHF  TO  DEGREES  CELSIUS. 

62  30! 


62^0  SUB  Tvsv ( V . T  ) 

6250  COM  /Co/  D ( 7 i . Aa ( 76.2) 

b260  Sun-0 

6270  FOR  I -0  TO  7 

6280  Sum*Sun*D ( 1 1 *V ‘  I 

6290  NEXT  I 

6300  T-Sun 

6310  SU8EN0 

6320! 

6330  !  CALIBRATES  T/C  READINGS. 


* , 3D . 3D, **  (H 
**,20.30,** 
**,2D. 3D.** 
** ,20. 3D, ** 
*“ .2D.30,  ** 

**.3D.DD,*:H 
**,3D.DD, *;H 
** .30.DD, * ;H 
** .3D.DD, * ;H 
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6340 ! 

63S0  DEF  FNTcor  r  <  T ,  I ) 

6360  COH  /Co/  0(7  I ,Aa<76.2> 
6370  Tc-Aa(I.O) 

6380  FOR  J-1  TO  2 
6330  Tc-Tc»rta( I, J ) *T  * J 
6400  NEXT  J 
64)0  RETURN  Tc 
6420  FNEND 


APPENDIX  B 
SAMPLE  CALCULATION 


Figure  (B.l)  shows  the  components  of  the  heat  balance  in  the  test  channel.  A 
sample  calculation  follows  that  details  the  methods  used  in  the  data  reduction  program 
to  estimate  the  value  of  the  energy  balance  components  and  then  calculate  the 
hydraulic  Reynold'  number  and  average  Nusselt  number  for  each  data  set.  A  copy  of 
an  actual  record  of  one  set  of  experimental  data  for  the  concave  test  section  precedes 
the  sample  calculations  together  with  a  list  of  the  parameters  necessary  to  make  the 
sample  calculation. 


Figure  B.l  Test  Section  Fnergy  Balance. 


SAMPLE  CALCULATION  DATA 


Month,  Day,  and  Time!  09:  16:  IS:  14:54 

THE  FOLLOWING  DATA  WERE  RECORDED: 

Heater  Vo  1  tape  tVM  *  91,75.5.  ( , 

Resistor  Volta ae  (Vr*  *  '.70 A  (VI 

•revision  Re- is  tor  <Rpr.  -  1.111729  >  Ohms. ) 


v,i.'  TEn  •  1  u  T  l  : 


T  ' hur.be  r  : 

1 

2 

u 

z 

6 

1 enp  ( v  ‘  : 

72.33 

t‘3.0" 

7  6 ’.00 

68 .01 

76".  30 

72.90 

* 'C  Nonce r: 

£• 

h 

!  0 

1  ! 

1  2 

Temp  ( C  1  1 

79.99 

72.13 

7. .77 

74.26 

77.03 

76.00 

T /C  Number : 

13 

1  4 

IS 

16 

1  ? 

IS 

T  emp  ( C i I 

74.32 

75.52 

75.70 

63.57 

62.45 

76.73 

T ''C  Number: 

13 

20 

21 

•>  '* 

23 

24 

T  emp  ( C i : 

73.13 

84.26 

8E.87 

69^58 

76.27 

76.07 

T /(  Number: 

TC 

26 

27 

28 

29 

30 

Temp  t  L  I  : 

77.53 

60.33 

47. 70 

79. 2S 

80.06 

83.47 

OUTER  INSULATION; 

T/C  Number: 

31  <  T i nsol  ) 

32  (Tinso2)  33 

l  T  inso3) 

T  emp  ( C l : 

82 

.38 

56.88 

40.01 

ORIFICE  TEMP 

(Tort)  • 

30. 25 

(C) 

INLET  TEMPERATURE: 

T/C  Number: 

35 

36 

37 

38 

T  emp  l C  1  : 

22.09 

22.07 

22.08 

22.07 

INNER  PLATE; 

T/C  Number: 

41 

42 

43 

44 

45 

46 

T  emp  ( C  I  : 

23.89 

23.69 

23.40 

23.34 

23.21 

23.49 

T/C  Number: 

47 

46 

49 

50 

51 

C 

T  emp  ( C  i  : 

23.28 

26.29 

23.56 

23.27 

24.03 

24~23 

T  /C  Number : 

53 

34 

99 

56 

57 

58 

T emp  ( C  >  : 

23.71 

22.82 

23.59 

23.63 

23.56 

24.02 

T/C  Number: 

S3 

60 

61 

62 

63 

64 

T  emp  i C  > : 

23.80 

24 . 49 

24.27 

24.87 

24.68 

25.87 

T/C  Number: 

65 

66 

67 

68 

63 

70 

T  emp  ( C  i  : 

25.56 

25 . 32 

2E .  1 6 

25.04 

24.95 

24.93 

INNER  INSULATION: 

T  /C  Number : 

71  ( T 

i  ns  i  1  * 

72  ( T i ns  i 

2)  73 

1  T i ns i 3 ) 

i emp  i C  )  : 

24 

.  4(i 

24.13 

23.70 

OUTLET  TEMPERATURE: 

T/C  Numbe r : 

74 

79 

75 

77 

T  emp  ( C  l  : 

30.66 

30.73 

30.65 

30.69 
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THE  FOLLOWING  TEMPERATURES  WERE  CALCULATED: 


Average  Outer  Wall  Temperature  (Two) 

Average  Inner  Wail  Temperature  (Twi) 

Average  Outlet  Temperature  (Tout) 

Average  Inlet  Temperature  (Tin) 

Channel  Inlet  and  Outlet  Temp  Diff  (Tdiff) 
Average  Bulk  Temperature  (Tblk) 

Mean  temperature  Difference  (Tdel) 

Outer  and  Inner  Wall  Temp  Difference  (Tudiff) 
Average  Local  Channel  Temp  at  Inlet  (Tcin) 
Average  Local  Channel  Temp  at  Middle  (Tcmid) 
Average  Local  Channel  Temp  at  Outlet  (Tcout) 
Avg  Local  Heated  Wall  Temp  at  Inlet  (Tumi 
Avg  Local  Heated  Wall  Temp  Near  Middle  (Tumid) 
Avg  Local  Heated  Wall  Temp  at  Outlet  (Tuout) 


7(5.46  (C) 
24. !6  (C) 
30.68  ( C ) 
22.08  < C ) 
3.61  (C) 
26.38  (C) 
SO. 08  (C) 

52.30  (C) 
22.44  (C) 

25.30  (C) 
30.32  (O 
73.62  (C) 
73.06  (O 
79.78  (C) 


THE  FOLLOWING  ORIFICE  DATA  WERE  ENTERED: 

Patrn(inHg)  DPM(inH20)  P!M(inH20>  RE1 

29.86  i 7 .60  16.90  1S000.0 

Dorf(m)  A(m"2)  Beta 

0.027318  S.861E-04  0.S277 


THE  FOLLOWING  DATA  WERE  CALCULATED: 

Orifice  Expansion  Factor  IY)  •  0.9859 

Orifice  Flow  Coefficient  (K)  -  0.6326 

Density  Based  on  Tort  < Rho >  •  1.1099  (kg/m'3) 

Viscosity  Based  on  Torf  (Mu)  *  18.61 6E- 06  (kg/m.s) 

Viscosity  Based  on  Tblk  (Mu)  *  13.434E-06  (kg/m.s) 

Therm  Cond  Based  on  Tblk  (Kair)  *  2.620E-02  (W/m.K) 

Patm  -  100848.13  (N/m‘2) 

Pdel  -  4379.58  (N/m'2) 

PI  -  96642.74  (N/m‘2) 

Mdot  •  0.0360  (kg/s) 

Repipe  *  4761 1 .3 

Red  *  7697.2 

Rehd  •  15018.9 

De  -  397.9 


op 

- 

350.411 

(Watts) 

Oai  r  • 

31 1 .991 

(Watts) 

Oiol 

m 

5.908 

<  Wa  t  ts ( 

Oiil  • 

.063 

(Watts) 

0 1  o2 

m 

3 .  '3u8 

(Watts) 

01 12  • 

.  1  On 

(Wattj) 

Olo 

m 

4.  30:3 

(  H  1  t  t  5  1 

Oli 

.  O0  1 

i  Wa  1 1  i 

it 

l  .  296 

(Matts i 

0oe  1  ■ 

30 .  wS3 

1 wa 1 1  s  i 

H  i  n 

m 

73. 7u 

( W/m ‘ 2C ) 

Nu  i  n 

•  3  7  .  t>  7 

Hm  i  d 

m 

84.33 

i  W/m  2<  ) 

Hum  i  d 

•  411.  J3 

Hou  t 

m 

31  .48 

• W/m ' 2C ) 

Nuou  t 

•  38.0  =3 

Havg 

30.4  7 

(  W/m ' 2C ) 

Nuavg 

•  38.06 

END  OF  RIJN 
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vh 

= 

91.755  V  Pwet  =  .5207  m 

V 

= 

7.704  V  r.  .297  m 

V 

= 

2.01729  SI 

p 

atm 

= 

29.86  in.  Hg 

AP 

= 

17.60  in.  H20 

Pi 

= 

16.90  in.  H20 

T,n 

- 

22.08  °C 

Tout 

= 

30.68  °C 

T 

wo 

= 

76.46  °C 

T  . 

wi 

= 

24.lt  °C 

Tinsil 

= 

24.40  °C 

T 

insi2 

= 

t^> 

OJ 

o 

n 

Tjnsi3 

= 

23.70  °C 

T 

insol 

= 

82.38  °C 

Tjnso2 

= 

56.88  °C 

TjnsoJ 

S 

40.01  °C 

4- 

k. 

O 

I- 

* 

30.25  °C 

Kair 

- 

.0262  W  m-°K 

CP 

- 

1007  J  Kg-°K 

^air 

4.6971  x  io-8  x  Taif  +  1.71947  x  l(f5  Kg 

Ki„s 

= 

.0380  W  m-°C 

AXin, 

.0127  m 

Ccu 

* 

0.12 

5.669  x  io'8  W  m2-°K4 

P 

- 

.5277 

Y 

- 

1.40 

- 

1.0  Kg-m  N-s2 

R 

- 

286.98  Nm  Kg-°K 

F 

WO-W| 

- 

1.0 

Apl 

- 

.07742  m2 

A 

c 

.0016  m2 

Ap'pe 

- 

.002 1 1  m2 

Dc 

* 

.00635  m 

Dpip« 

- 

.<>518  m 

Dorf 

- 

.027318  m 
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2.  TEMPERATURE  CALCULATIONS 


a.  Average  Temperature  of  the  Outer  (Inner)  Wall 

_  _  V  Outer  (Inner)  Plate  Temperatures 

Two(Twi)  "  r - r - — - 77 - - - — - - -  =  76.46°C(24.16°C) 

Number  of  Outer  (Inner)  Plate  Thermocouples 

b.  Average  Channel  Outlet  (Inlet)  Temperature 

-r  /-T-  v  X  Outlet  (Inlet)  Temperatures 
out  in  "  Number  of  Outlet  (Inlet)  Thermocouples  =  3a68°C  (22-08°C) 

c.  Channel  Inlet  to  Outlet  Temperature  Difference 

TdifT  =  Tout  -  Tin  -  30-68°C  -  22-08°C  =  8.6°C 


d.  Bulk  Temperature  of  the  Flow  Through  the  Channel 


Tb.k  = 


T.  +  T 


out  _ 


22.08°C  +  30.6S’C 


=  26.38°C 


e.  Average  Temperature  Difference  of  the  Heated  Wall  and  Fluid  Bulk  Temperature 

AT  -  Two  -  Tb]k  =  76.46°C  -  26.38°C  =  50.08°C 

f.  Temperature  Difference  Between  the  Outer  and  Inner  Walls 

Twdiff  *  Two  *  Twi  =  76  46’C  -  24.16’C  =  52.30‘C 


3.  POWER  CALCULATIONS 


a.  Power  supplied  by  the  Heated  Plate 

V.xV  91.755V  x  7.704V 
Q>  ■  -  V*  -  -  2.0172911  ■  550  41  Wam 


b.  Power  Loss  Through  the  Outer  Wall  Insulation 
Q,»  ‘  «  2«Q.oi  + 

“  (l/2K(^^W)timA|)| 

ins 

T  -  T 

4.  t  inso2  inso3  \i.  a  i  _ 

AX  )kinsApli 

ins 

,  ,  (0.038  W  m°C) 

(1  2)(. 07742m2) - — — - ((82.38  -  56.88(°C))  +  (56.88  -  40.01(°C))]  =  4.908W 

.0127m 

c.  Power  Loss  Through  the  Inner  Wall  Insulation 

Q,i  =  (i'2KQ m  +  Qli2) 

=  (1  2)((— S"T|nSi2  )RinsApi 
ins 


T  -  T 

insi2  insi3 


■>kinsApll  = 


,  (0.038  W  m’C) 

(1  2)(. 07742m*) - — — - [(24.40  -  24.13(°C))  +  24.13  -  23.70(°C))]  =  .0S1W 

.0127m 

d.  Heat  Loss  by  Radiation 

/.  Radiation  Resistance 
1  -  e  1 

R  =  2  x  ( - £U_)  +  -  = 

r  AP.  X  £cu  Ap.  *  *wo  wi 

12  12  , 

=  - - - 1)  =  - ,( - 1)  =  202.36  m2 

Ap,  Ccu  0.07742m2  0.12 

2.  Heat  Radiated 
cs  -g(Two4-V) 


(5.669  x  io  8 >((76.46  + 273 )4  -  (24.16+  273)4] 


202.36 


=  1.996  W 


r*  T  t, 


4.  MASS  FLOW  RATE  CALCULATIONS 


4! 


$ 

!■# 


-*<! 

V*' 


,v 

.»!J 


a.  Pressure  Conversions 

Paim  =  29  865  in  h8  x  3376.8  =  100848.13  \  m2 

AP  =  17.60  in  H20  x  248.84  =  4379.58  N  m2 

P,  =  (16.90  in  H,0  x  248.84)-  100S48.13  *  96642.74  N  m2 

b.  Density  of  Air 

P.  (96642.74) 

Pair  RTorf  ~  (286.98)  (30.25  +  273)  *  1-1099  Kg  m 

c.  Expansion  Factor 


Y  =  1  -  (0.41  +  0.35p4) 


=  1  -  (0.41  +  0.35(0.5277)4] 


d.  Mass  Flow  Rate 

m  =  YKA  J1TS..A'P 


4379.58 

(1.401(96642.74) 


=  0.9859 


m  =  YKA  V2gcpairAP  =  (.9859)K(.0005S61m‘) 

x  ^2(1  Kg-m.Ns2X  1.1099  Kg  nv)(4379.58  N  m2>  -  .05669 K  Kg  s 
Iterate  to  find  m: 

Assume  a  hydraulic  Reynolds  number  *  15000 
D.  A  .05 1 8  x. 00 16 

pipe  DhdApipe  hd  .01229  x. 002 II  hd 

-»  Repipe  ~  15000(3.1961)  =  47942 

Interpolate  in  the  appropriate  table  in  Reference  15  for  K.  the  flow  coefficient: 

K  =  .6329 

(.Vo/e  that  the  c<  nputer  program  used  to  reduce  the  data  used  an  appre.ximai.ne 
polynomial  to  Jind  e  correct  value  of  the  jlow  coefficient,  K.) 

-*  m  =  (. 05669 )(. 6329)  =  .0359  Kg  s 
Solve  for  the  improved  value  of  Re 

r  pipe 

_  m  D.  .  (.0359X.O51S) 

Re  .  =  - BJCS_  =  - - -  =  4",|S 

ApipeMair  (.W21 1  K  18.616  x  1  ( ff' ) 

Iterate  until  successive  values  of  m  differ  h\  le<<;  than  <*>1  m  =  » 


5.  REYNOLDS  NUMBER  CALCULATIONS 


a.  Re  Based  on  Pipe  Diameter 

(.0360X.0518) 


R.,„  -  = 


P*Pe  A  u  (.00211X18.616  x  10"6) 

pipe'air  v  7 

b.  Re  Based  on  Channel  Height 

(.0360X. 00635) 


=  47475 


m  D 

Re, - = 


d  Ac^air  (.0016X18.434  X  |0'6) 


=  7751 


c.  Re  Based  on  Hydraulic  Diameter  of  the  Channel 


Rehd  = 


m  D 


hd  _ 


A  It 

c^air 


(.0360X  01229) 
(.0016X18.434  x  lo~b) 


=  15000.8 


6.  HEAT  TRANSFERRED  TO  THE  WORKING  FLUID 


Qa,r  =  n’KP  (Tou,  ’  Tin)  =  (  °360  K?  sK  1007  J  Kg-eC)( 30.86  -  22.08 >X  =  312  W 


AVERAGE  HEAT  TRANSFER  COEFFICIENT  CALCULATION 


312  W 


Ap,AT  (. 07742m2  )(50.085C) 


=  80.5  W  m2-9C 


8.  AVERAGE  NUSSELT  NUMBER  CALCULATION 


_  h  Dhd  (80.5  W  m‘-9C)(. 01 229m-) 


N 


“hd 


K. 


air 


.0262  W  m-'C 


=  37.8 


9.  DEAN  NUMBER  CALCULATION 
Re 


De  = 


-  V 


/D  2  7T5l  oo<> 3 5~ 


V 


=  400.7 


"TVW»  ■ 


APPENDIX  C 

EXPERIMENTAL  UNCERTAINTY 

The  experimental  uncertainty  in  calculated  parameters  was  deternuned  in 
accordance  with  the  method  described  by  Kline  and  McChntock  [ Ref .  2?) 

The  qualitative  results  shown  here  correspond  to  the  data  set  and  calculation  ol 
Appendix  B.  There  was  a  significant  decrease  in  the  experimental  uncertatntv  of  the 
average  Nusselt  number  for  Reynolds  numbers  above  I'.OOO  in  this  experimentation 
compared  to  that  of  Galyo  [Ref.  14:  pp.  104-105]  because  of  the  changes  made  to  'he 
test  channel  that  promoted  mixing  of  the  working  fluid  at  the  exit  ol  the  channel. 

The  equations  used  to  determine  the  experimental  uncertainty  of  calculated 
parameters  are  shown  below. 

1.  DENSITY  OF  THE  WORKING  FLUID 

8p  6R  ,  8P.  ,  6T  ,  ,  n , 

P  R  Pl  Torf 


2.  MASS  FLOW  RATE  OF  THE  WORKING  FLUID 

6m  SY  ,  5K  ,  6A  , 

— -11'— >  +  <  — >  +  1 

+  (I  4)((  5£ai£ —  )2  + 

p  AP 

^air 


3.  HEAT  TRANSFERRED  TO  THE  WORKING  FLUID 

5Q  .  5m  ,  6c  , 

=K  —  )2  +  (  — )2 

Q,ir  m  C 


6<ToU.  '  Tin) 

<Tout '  Tin) 


2  i  0.5 
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4.  AV  ERAGE  heat  transfer  coefficient 


5L 


l( 


*Q. 


^nr 


( >2  +  ( )2  J° 5 

Ap.  at 


average  nusselt  number 


4Nil.  , 

p»d-  -  l( - r  + 


6D 


D 


r  *  ( 


,0! 


hd 


air 


HYDRAULIC  REYNOLDS  NUMBER 


$Re 

~rT 


dm  , 

l<  —  r 

m 


6D 


D 


M—  ,2 


hd 


air 


6AC 

A 


>2  l05 


Several  parameters  used  in  the  calculations  described  herein  had  uncertainties 
that  were  constant  over  the  entire  range  of  testing.  These  quantities  and  their  relative 
uncertainties  follow 


Quantity  <  B ) 


Uncertainty  (  ±6B  B) 


A 

A 

A 


c 

pi 


p 

D 


D, 

K 

K 


hd 


air 


atm 


.0020 

(>408 

.0006 

0020 

TMOO 

.0391 

.0030 

(«)38 

.0010 

(8X)3 

<8180 

.0020 

<8)54 
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The  approximating  polynomials  used  to  resolve  thermocouple  voltages  into 
temperatures  matched  actual  temperatures  to  within  less  than  ±0  1C  over  a  range  of 
temperatures  from  204C  to  85  C;  however,  in  some  cases  this  temperature  range  was 
exceeded.  A  conservative  estimate  of  the  uncertainty  of  temperatures  measured  bv 
thermocouples  was  set  at  ±0  1C.  the  accuracy  of  the  thermocouples  specified  by  the 


manufacturer  of  the  thermocouple  wire, 
calculations: 

The  following  uncertainties  were  used  in 

Quantity  (T) 

l  ncertaintv  (  ±  T) 

V 

o.l  C 

T,n 

o.l  C 

Torf 

o.l  C 

o.l  C 

u 

O.l  C 

Tout  *  Tm 

0.2  C 

AT 

0.2  C 

For  the  data  taken  at  a  hydraulic  Revnolds  number  of  15000  in  the  concave  test 
section  the  following  additional  uncertainties  applied: 

Quantity  <  P»  l  ncertainty  (  ±  Pi 


P|  Dll' 

AP  o  |ii 


The  results  of  the  uncertainty  calculations  at  a  Reynolds  number  of  15000  in  the 
concave  test  section  were: 


Quantity  (B) 


Uncertainty  ( ±6B  B) 


p 

.00445 

m 

.00544 

Qair 

.02397 

h 

.02431 

.04620 

Rehd 

.05703 

The  major  contribution  to  experimental  uncertainty  in  both  the  average  Nusselt 
number  and  hydraulic  Reynolds  number  was  the  uncertainty  in  cross  sectional  area  of 
the  channel.  The  hydraulic  diameter  of  the  channel  was  a  function  of  channel  cross 
sectional  area,  resulting  in  a  relatively  large  uncertainty  in  the  hydraulic  diameter. 
These  two  parameters  accounted  for  approximately  fourty  percent  of  the  uncertainty  in 
the  average  Nusselt  number  calculation,  and  about  eighty  percent  of  the  uncertainty  in 
the  calculated  hydraulic  Reynolds  number. 
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